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Preface 


The  fifth  Forum  on  Space  Structures,  like  its  four  predecessors,  sought  to 
promote  intense  dialog  on  contemporary  issues  in  the  interaction  of  controls 
and  dynamics  for  large  flexible  spacecraft.  To  this  end  the  meeting  followed 
the  forum  tradition  of  open  discussions  interspersed  by  short  (10  minute) 
presentations . 

The  Forum  on  Space  Structures  originated  from  a  need,  perceived  by  a  group  of 
Air  Force  supported  research  scientists  and  engineers,  for  improved  communica¬ 
tion,  interaction,  and  col  1 aoorat i on  between  the  structural  dynamics  and 
controls  disciplines.  They  considered  this  essential  if  the  well  recognized 
potential  benefits  of  active  control  of  dynamic  response  were  to  be  realized. 
Consequently,  a  group  of  these  individuals  gathered  in  an  informal  meeting  at 
MIT  in  1983.  The  free  exchange  of  ideas  that  material i zed  fron  the  open 
discussion  format  adopted,  was  hailed  by  all  participants  as  worthy  of 
repetition  at  least  annually  for  the  next  several  years.  The  subsequent 
neetings  have  continued  to  foster  increased  communication  and  interaction 
between  the  controls  and  dynamics  communities.  The  fifth  Forum  was  no  less 
successful  of  Us  predecessors. 

This  document  is  an  attempt  to  record  the  essence  of  the  discussions  for  the 
benefit  of  other  interested  readers  in  the  hope  of  extending  whatever  accomp¬ 
lishments  were  attained  beyond  the  specific  space  and  time  confines  of  the 
meeting.  Abstracts  of  the  brief  presentations  have  been  reproduced  with  only 
minor  editing.  In  addition,  summary  accounts  of  the  discussions  acconpanying 
and  "hollowing  the  presentations,  as  recalled  by  several  participants,  have 
been  compiled  and  included  herein.  An  executive  summary  has  also  been 
provided  that  summarizes  the  major  highlights  and  achievements  of  the  fifth 
Forum. 

On  behalf  of  the  Air  Force,  we  wish  to  express  our  appreciation  to  all  who 
participated  in  the  discussions,  and  in  particular  'Jr  Robert  Kosut  who  took 
on  the  very  difficult  task  of  organizing  and  chairing  the  meeting,  and  did 
them  both  so  well.  We  are  also  grateful  to  Professor  Walter  Pi  1  key  who  also 
accepted  the  burden  of  overseeing  the  production  of  this  document.  Lastly, 
we  commend  those  participants  who  volunteered  to  and  did  taxe  notes  of  the 
proceedings,  and  subsequently  composed  the  accounts  contained  herein. 
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Synopsis 


The  issues  of  discussion  at  the  Fifth  Forum  pr  vided  clear  indications  of  progress 
made  in  the  Forum  objectives  since  the  inception  of  the  meetings.  Technical  progress 
was  evidenced  by  the  transition  of  the  focus  from  the  topics  that  used  to  dominate  the 
discussions  to  those  that  were  then  considered  too  esoteric  and  too  far  out  to  be 
meaningfully  addressed  by  the  audience. 

In  the  modeling  area,  the  primary  concerns  used  to  be  on  issues  related  to  the 
synthesis  of  controllers  for  complex  flexible  structures.  They  included  basic  modeling 
approaches  -  distributed  versus  discrete  parameterization,  accommodation  of  parameter 
uncertainties  and  modeling  errors,  model  reduction  techniques  and  criteria,  and  system 
identification  for  model  updating  in  both  synthesis  and  adaptation.  In  contrast,  the 
modeling  issues  of  the  Fifth  Forum,  with  but  one  exception,  dealt  principally  with 
performance  representation  and  prediction,  and  with  performance  implications  of 
actuator  dynamics.  The  one  exception  was  the  Integrated  Structural  Analysis  and 
Control  (ISAAC)  presentation  which  highlighted  progress  to  date. 

The  significance  of  orbital  effects  such  as  gravity  gradient  and  dynamic  coupling 
between  rigid  and  vibratory  motions  in  the  dynamics  and  stability  of  a  class  of 
spacecraft  systems  was  the  issue  in  Peter  Bainum’s  presentation  and  accompanying 
discussions.  These  effects  were  also  the  motivating  considerations  in  the  modeling 
approach  outlined  by  Lcn  Meirovitch.  Ed  Haug’s  presentation  also  dealt  with  this 
general  theme  but  from  a  simulation  and  computational  algorithm  standpoint.  Jim 
Williams  led  a  very  lively  discussion  on  travelling  waves  and  their  Telcvance  to  orbital 
performance.  This  is  perhaps  one  of  the  perennial  topics  of  interest  at  these  meetings. 
But  in  spite  of  its  popularity  only  a  few  attempts  have  been  made  to  incorporate  wave 
mechanics  approaches  in  the  modeling  of  large  spacecraft  dynamics.  Dan  Inman  and 
Walter  Pilkcy  both  addressed  the  issue  of  actuator  dynamics  and  their  impact  on  system 
performance  (damping).  Finally  the  issue  of  in-flight  adaptation,  well  recognized  for 
control  systems,  was  broadened  by  Michail  Zak  to  include  structural  parameters,  c.g. 
stiffness. 

In  the  controls  area,  the  shift  in  focus  was  more  subtle,  with  controller  quality  issues 
such  as  robustness,  stability,  and  implementability  dominating  over  the  purely 
algorithmic  issues.  The  wave  motion  control  issue  raised  by  David  Miller  represented  a 
significant  change  from  previous  discussions  in  which  the  emphasis  was  clearly  on 
modeling  considerations.  Robert  Kosut’s  discussion  of  Adaptive  Controls  provided  a 
brief  review  of  the  state  of  development  in  modeling,  parameter  identification, 
robustness,  stability  and  convergence.  Stability  analysis  is  currently  feasible  for  slow 
adaptation,  and  extensions  to  faster  adaptation  rates  represents  the  current  thrust  of 
research.  In  a  slight  departure  from  his  prepared  abstracts,  Dave  Hyland  discussed  issues 
from  thcorv/practicc  considerations,  including  difficulty  of  modeling  joint  nonlincaritics 
and  local  resonance  conditions,  both  of  which  played  a  major  role  on  the  COFS-I  project. 
He  also  described  some  experimental  development  programs  at  Harris  including  an 
analytical/cxperimcntal  program  being  initiated  under  SDIO/IST  funding.  The  object  is 
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to  develop  and  validate  a  hybrid  design  qualification  methodology.  His  two  abstracts 
dealt  with  the  determination  of  performance  bounds  caused  by  modeling  uncertainties. 
Decentralized  control  algorithm  developments  were  discussed  by  Umit  Ozguner,  David 
Young  and  Juri  Medanic.  The  first  also  included  a  discussion  of  experimental  programs 
conducted  or  planned  to  validate  the  analytical  developments.  David  Young’s 
presentation  was  particularly  interesting  because  it  approached  the  decentralized  control 
problem  from  the  context  of  physical  components  of  the  system  with  local  control 
requirements  and  arrived  naturally  at  a  hierarchical  (interlocking)  configuration  for  the 
global  system  controller.  It  provided  reinforcement  to  the  effort  at  MIT  by  Ed  Crawley 
in  which  the  hierarchical  concept  was  assumed  by  physical  intuition.  John  Junkins 
concluded  the  formal  presentations  with  a  discussion  of  his  multi-criteria  optimization 
scheme  for  robust,  minimum  control  effort  designs  based  on  eigenvalue  assignment. 

The  concluding  general  discussion  touched  on  many  of  the  issues  raised  during  the 
meeting  but  focussed  principally  on  the  problems  of  validation  of  theory,  and  on 
transition  to  application.  On  the  latter  it  was  felt  that  transition  is  lagging  the  advanced 
technology  developments  by  about  30  years.  This  was  considered  unsatisfactory  for  the 
structurc/controls  interaction  disciplines.  Some  counter  examples  were  however  alluded 
to  in  the  helicopter  and  fighter  aircraft  areas.  On  the  theory  validation  issue  there  was  a 
strong  consensus  on  the  need  for  flight  experimentation  to  provide  the  feedback 
necessary  for  fine-tuning  of  the  theories  prior  to  transition.  It  was  observed  without 
further  discussion  that  current  funding  trend  seems  to  favor  experimental  over 
theoretical  research  projects. 

Progress  in  the  nontechnical  objective  of  the  Forum  -  promotion  of  improved 
communication  between  the  Controls  and  Structures  communities  -  was  also  evident  in 
the  degree  of  interaction  that  occurred.  No  longer  were  the  two  communities  directing 
questions  at  each  other  from  opposite  camps  or  pointing  accusative  fingers  at  each  other 
on  specific  issues,  but  all  participated  actively  in  addressing  the  various  topics  each 
contributing  to  the  discussion  from  the  vantage  point  of  their  own  expertise. 
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Modeling  of  Flexible  Spacecraft 
Accounf'ng  for  Orbital  Effects 
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INTRODUCTION 

Mnnv  current  investigations  oi  tlio  shape  iml  ori-m!  ion  control  oi 
proposed  I  lex i bio  orbit  in*  l a  rue  space  structures  ( I.SS  •  To  not  i m  orpo- 
rate  the  effects  of  the  uravit  v— gradient  and  orbital  dynamic  coupl i  ng 
into  the  plant  mode  1  s  .  ’i'll  i  s  means  that  for  the  corresponding  l  ino  ir  i  zed 
unforced,  open-loop  systems  the  poles  of  the  rigid  rotational  modes  are 
at  the  origin.  The  manner  in  which  the  orbital  1 v  induced  coupling  ef¬ 
fects,  due  to  gravity-gradient  ind  -gyroscopic  effects,  are  introduced  is 
clearly  indicated  in  the  continuum  lonnulution  o!  Saittini  1  lor  pre¬ 
dicting  the  motion  ot  a  general  orbiting  i lexible  bodv  in  orbit.  These 
coupling  terms  reflect  both  coupling  between  the  rigid  and  1 lexible  mo¬ 
tions  ind  also  inrra-mod.il  coupling  attests.  Plastic  deiormnt ions  ire 
considered  small  as  compared  with  charac ter  i  si  ic  bodv  dimensions.  liqua¬ 
tions  are  developed  tor  both  the  rigid  and  elastic  (generic)  motions, 
based  on  m  i  priori  knowledge  ot  i he  frequencies  and  shape  functions  ot 
all  modes  included  within  the  truncated  system  model.  The  orbital  Ly 
induced  coupling  terms  are  seen  to  depend  on  volume  integrals  whose  in¬ 
tegrands  i  re  : unctions  ot  r  lie  various  components  ol  the  dilferenr  modal 
shape  1  unctions  together  with  the  coordinates  ol  the  dilierenti.il  mass 
e  1  erne  nf.s  1  ,  . 

NUMERICAL  STUDIES 

This  paper  summarizes  two  recent  numerical  studies  aimed  at  attempt¬ 
ing  to  evaluate  the  el  leer,  ot  i nc 1 nd i ng/ omi l t i ng  the  orbit  il  .oupling 
terms  in  the  modelling  ol  I.SS  systems.  The  first  is  based  on  a  mode  1  •'! 
an  orbiting  shallow  spherical  shell  svston  whieh  com  lins  >  hvnrid  sys¬ 
tem  lor  controlling  both  the  shape  ind  the  oriental  i  on  (Fig.  1  i.  The 
passive  sp r i ng- 1 oaded  dumbbell  type  damper  att  idled  at  the  shell's  ipex 
is  ilso  used  i  o  provide  the  l.ivorihle  moment  ot  inert  i  .  ri;  it'  required 
f  or  .;  r  i  v  l  ;  v-gr  oi  i  en-  .tab  i  1  i  cat  i  on  .  don-frol  laws  ir,  ve  i  epe.l  !  >r  tie- 
six  shell  mount  ed  point  ictuitors  based  <>n  pole  cluster  i  ng  t  ecliniqiirs. 

A  studv  was  mute  to  determine  the  el  tecta  o!  omitting  the-  coupling  t  erms 
in  the  development  ot  control  1  iws.  lont  rol  laws  based  on  pole  pi  icemont 
were  t  irsf  developed  tiased  on  the  shell  model  which  d  i  d  not  include  tile 
grav  i  t  v-gr  id  i  ent  and  orbit  il  coupling  terms  in  the  pi  mt  model.  fhe  ioii- 
trol  laws  thus  developed  were  then  inserted  inti)  the  previously'  aevel- 
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oped  models  which  contain  both  iirst  order  gravity-gradient  and  orbital 
dynamic  coupling  terms.  Typical  results  are  illustrated  in  Fig.  2  and 
Ref .  [ 2].  In  all  cases  studied,  there  is  a  general  tendency  tor  some  of 
the  poles  ot  the  rigid  modes  to  shiit  towards  the  imaginary  axis  when 
the  gravity-gradient  and  orbital  effects  are  superimposed  into  t’ae  plant  . 
For  the  less  robust  systems  (Fig.  2)  instability  may  result.  In  gemril, 
there  is  no  noticeable  shift  in  the  poles  corresponding  to  the  flexible 
modes.  As  expected,  for  the  more  robust  systems,  the  relative  effect  ot 
this  shift  is  less  apparent  and  stability  may  still  be  maintained,  but 
at  the  expense  of  greater  control  force  eliort  ;2]. 

In  the  second  study  a  computational  algorithm  was  developed  to  eval¬ 
uate  the  coefficients  of  the  various  coupling  terms  in  the  equations  ot 
motion  as  applied  to  the  finite  element  model  ( FFM  )  ot  the  122m  diameter 
Hoop/Column  system  t  ):.  For  the  discretized  FFM  the  expressions  for  the 
volume  integrals  in  the  coupling  coelficionts  wore  replaced  by  summations 
of  similar  terms  taken  over  the  total  number  of  discrete  masses  for  which 
eigenvector  (mode  shape)  information  was  available  from  the  FEM  output. 
The  details  concerning  these  calculations  are  provided  in  Ret.  3'. 

Table  1  shows  the  comparison  between  the  magnitudes  ot  the  components 
ot  the  main  coel f icients  in  Euler’s  rigid  rotational  equations,  R,  ind 
the  components  of  the  largest  coupling  coe  I  1"  i  c  i .  n  t  s  ,  Q  .  for  the  seven 
flexible  modes  contained  in  the  truncated  FEM  model.  Two  different,  val¬ 
ues  of  modal  amplitude  factors,  An,  of  1m,  ind  1mm.  respectively,  were 
assumed  lor  all  the  flexible  modes.  It  appears  that  when  the  system  is 
operating  within  the  mission  specifications  (deflections  of  the  order  ot 
mm),  the  finite  element  assumptions  neglecting  the  orbital  coupling  are 
valid;  when  the  deflections  are  of  the  order  of  meters  the  coupling  be¬ 
tween  the  flexible  and  rigid  modes  should  he  incorporated  into  the  equa¬ 
tions  of  mo t i on . 

In  Table  2.  a  comparison  of  pertinent  terms  in  the  generic  modal 
equ  it  ions  is  given.  Rn  An  *  w“A„  is  compared  with,  coupling  terms  repre¬ 
sented  by  !’n.  The  comparison,  after  normalization  of  the  different  terms 
in  the  equations,  shows  that  the  time  dependent  amplitude  ot  t  lie  modes 
can  be  approximated  as  an  harmonic  oscillator  at  least  up  to  the  point 
where  An-lm  tor  all  n. 
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CONCLUDING  COMMENTS 

Even  though  tli''  numeric  1 1  in.il  vs  is  his  b»  •  <.•  n  ippLied  to  two  spec! t 
mo  d  «._•  Is  oi  proposed  orbit  inc  I.SS  svstiOb .  r  Me  results  indicate  that  one 
should  not.  .  a  priori,  neple.t  orbit  illv  induced  coup  l  inc  terms  in  the 
malvsis  mu  desi.cn  oi  Open  and  closed- loop  hSS  systems  in  general. 
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In  response  to  comments  raised  concerning  the  two  examples 
presented  and  the  importance  m-iai  r  tiiese  studies  the  followin 
additional  points  should  be  empha s  i  zed  : 

1)  For  both  sets  of  numerical  examples  presented 

it  appears  that  the  influence  of  orbitally  induced 
coupling  on  the  rigid  body  motions  can  be  more 
important  than  the  similarlv  induced  effects  on  the 
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disturbance  torques  on  the  vibrating  and  *  he  i  i  I  » 
deflected  system.  ('n  1  v  then  can  final  conclusions 
be  made  cence  rn i ng  RMS  pointing  and  shape  re. pii  fo¬ 
ments.  1  iu'  studies  cited  in  Refs,  2  and  3,  directed 
and  supported  by  NASA  l.anglev,  were  intended  to  gain 
ins  i  gilt  into  a  complex  problem  and  not  to  present 
a  final  design  of  t  tie  control  system  for  an  actual 
LSS  flight  svstem,  which  will  probably  depend  >n 
future  development-  in  i.SS  sensor  technology. 

4)  Regarding  slew  maneuvers,  as  long  as  the  elastic 
deformations  can  he  maintained  small  as  compared 
with  elia  rac  t  e  r  i  s  t  i  e  body  mule  formed  dimensions 
(such  as  length  or  diametorl  the  techniques  of  Ref. 

3  could  also  be  used  to  evaluate  the  coupling  terms 
during  or  immediately  following  a  slew  maneuver. 

For  violent  slews  where  large  elastic  deformations 
may  result  the  nature  of  the  coupling  terms  would 
change  from  those  of  Rot's.  1  and  3  and  a  different 
formulation  to  represent  the  dynamics  would  be 
required.  In  connection  with  slewing,  it  should  be 
noted  that  several  proposed  I.SS  systems,  such  as 
the  Wrap  Rib  Antenna  and  t  iie  sCPl.H ,  involve  asym¬ 
metrical,  of! set  sv<: ems .  desired  orientations  after 
line-  f  •  rare-  :\  e  r  -  "a.  :-.<t  a  1  ways  correspond  to 

gravitational  • '  •  r  i  i  !  l  i  • :  :  ;;  .  ui  d  it  ions,  so  that  oi;!:vi 
requirements  fo i  slewing  mav  vary  great  1 »  depending 
on  the  sense  or  direction  of  the  slew. 

5)  One  member  of  the  audience  raised  questions  about 
the  need  for  the  types  of  studies  described  in  Ref. 

2  and  3  since  lie  thought  that  the  effects  of  gravity 
induced  effects  on  the  closed  loon  dynamics  would  be 
small.  The  author  agrees  that  with  a  robust  enough , 
properly  designed,  and.  functioning  (no  fault!  control 
and  sensor  svstem  those  effects  can  he  properly  com¬ 
pensated  for  (big.  <-  of  Ref.  2).  In  the  event  o; 
failure  of  a  part  •  or  .ill  of  the  control  sens<  r  s  y  s  - 
tern  the  open  -  loop  lesult.  •  f  Re  f  .  3  might  also  ;•  •  •  v  ;  J 
usetul.  In  conclusion  the  author  would  like  te  t  hails 
the  mane  participants  f.  r  *  ho  i  r  interest  ..ad  help!:;!, 
c  on  st  met  i  v  e  c  ■'".me  n  t  •-  . 
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INTRODUCTION 

The  first  step  in  the  mathematical  analysis  of  any  physical  system  is  the  selection  of  a 
mathematical  model  to  represent  the  system.  This  selection  is  extremely  important,  since  it 
not  only  determines  in  advance  the  scope  of  possible  results  of  the  analysis,  but  also  heavily 
influences  the  design  of  auxiliary  devices  and  systems,  such  as  control  systems.  In  the 
dynamic  analysis  of  large  space  structures,  mathematical  models  have  consisted  primarily 
of  a  set  of  vibration  modes  [1,2].  The  popularity  of  the  modal  vibration  model  of  large  space 
structures  is  due  partly  to  the  well-developed  analytical  techniques  which  can  be  applied 
to  such  a  model,  and  partly  to  the  success  with  which  modal  vibration  models  have  been 
used  to  describe  structures  on  earth.  However,  in  view  of  the  unprecedented  size  of  large 
space  structures  and  their  potential  technological  importance,  it  is  worthwhile  to  examine 
the  limitations  of  modal  vibration  models  and  to  consider  the  usefulness  of  other  models.  In 
particular,  this  paper  is  concerned  with  models  which  view  large  space  structures  as  media 
in  which  wave  propagation  provides  an  accurate  dynamic  description. 

WAVES  AND  MODES 

The  relationship  between  modal  analysis  and  wave  propagation  analysis  may  be  illus¬ 
trated  by  the  Venn  diagram  in  Fig.  1.  The  relationship  is  analagous  to  the  relationship 
between  linear  elasticity  and  non¬ 
linear  elastic-plastic  continuum 
theory.  For  linear  problems,  wave 
propagation  analysis  and  modal 
analysis  are  mathematically  equiv¬ 
alent,  but  the  two  approaches  of¬ 
ten  give  enlightening  complemen¬ 
tary  perspectives.  Wave  analysis 
is  applicable  to  many  problems  in 
which  the  concept  of  modal  super¬ 
position  cannot  be  applied,  such 
as  the  problem  of  dynamic  failure 
and  other  nonlinear  problems.  In 
many  linear  problems,  such  as  the 
linear  analysis  of  joint  dynamics  or  the  development  of  linear  control  systems,  wave  propa¬ 
gation  analysis  may  have  computational  or  conceptual  advantages. 


Fig.  1  Relationship  between  modal  analysis  and  wave 
analysis. 
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Linear  Analysis 

In  linear  analysis,  when  wave  propagation  analysis  and  modal  analysis  are  theoretically 
equivalent,  the  preference  between  a  wave  analysis  and  a  modal  analysis  is  assessed  in  ac¬ 
cordance  with  the  ratio  between  the  characteristic  structural  length  L  and  a  characteristic 
length  associated  with  a  typical  disturbance  or  loading  function.  The  characteristic  length 
associated  with  the  disturbance  function  may  be  the  wavelength  A  introduced  into  the  struc¬ 
ture  by  a  point  disturbance  source,  the  spatial  extent  over  which  a  distributed  disturbance 
source  is  applied,  or  the  spatial  extent  over  which  a  disturbance  from  a  point  source  is 
effectively  attenuated  as  it  propagates  through  the  structure.  For  the  case  of  a  point  source 
disturbance  in  a  nondissipative  structure,  three  characteristic  regimes  of  dynamic  behavior 
may  be  defined.  When  A  >>  L,  a  quasistatic  or  rigid  body  dynamic  analysis  is  likely  to 
be  adequate.  For  A  ss  L,  the  structure  begins  to  deform  dynamically,  and  a  modal  analysis 
consisting  of  globally  defined  modes  is  likely  to  be  preferable.  For  A  <<  L,  the  globally  de¬ 
fined  modes  cannot  easily  represent  the  short  wavelength  disturbances,  and  a  wave  analysis 
is  likely  to  become  preferable. 

Structural  Details 

For  disturbances  with  short  wavelengths  (A  <  <  L),  structural  details  such  as  individual 
members  and  joints  become  important  in  predicting  dynamic  response.  An  advantage  of 
wave  propagation  analysis  is  that  it  can  deal  with  local  models  of  structural  subsystems.  By 
considering  transmission  of  structural  disturbances  through  a  subsystem,  a  wave  analysis 
can  be  based  only  on  a  local  model  of  the  subsystem  and  the  adjacent  portions  of  the 
structure  to  which  it  is  attached.  A  modal  analysis,  on  the  other  hand,  must  be  based  on 
a  model  of  the  entire  structure,  and  requires  a  large  number  of  accurately  modeled  modes 
in  order  to  include  the  significance  of  structural  details. 

Time  Delays 

Recent  results  in  the  analysis  of  large  tetrahedral  lattice  structures  show  that  the  lime 
delay  between  an  input  at  one  location  in  the  structure  and  the  response  at  another  location 
may  be  on  the  order  of  tens  of  seconds  [3j .  Such  time  delays,  which  may  have  important 
implications  for  control  system  design,  are  an  essential  part  of  a  wave  propagation  analysis. 

Dynamic  Failure 

Consideration  of  structural  failure  is  an  important  aspect  of  the  reliable  implementation 
of  large  space  structures.  In  a  dynamic  failure  problem,  the  geometry  and/or  the  material 
properties  of  a  structure  vary  with  time.  Thus,  a  normal  modal  analysis  cannot  be  applied. 
A  wave  propagation  analysis  can  be  used  to  predict  the  conditions  under  which  dynamic 
failure  will  occur  and  to  suggest  some  means  of  achieving  failure  arrest  1,  5,.  The  dynamic 
failure  problem  is  an  example  of  a  problem  which  lias  no  modal  desc’ription,  but  which  can 
be  analyzed  using  wave  propagation  concepts. 

CONCLUSION 

The  final  judgement  of  any  mathematical  model  must  be  based  to  a  large  extent  on  the 
ability  of  the  model  to  predict  those  phenomena  of  interest.  A  good  model  is  a  model  which 
may  be  used  to  design  and  operate  a  useful  physical  system.  Since  large  space  structures 
are  envisioned  to  be  complex  structures  which  are  to  be  used  for  a  wide  variety  of  purposes, 
it  is  unlikely  that  any  single  model  (a  modal  model  or  a  wave  model)  will  be  sufficient 
for  all  aspects  of  large  space  structure  applications.  Thus,  it  seems  unwise  at  this  point 
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to  restrict  attention  to  any  single  viewpoint.  When  large  space  structures  come  closer  to 
actual  deployment,  operational  reality  alone  will  determine  the  relative  merit  of  any  model 
which  is  considered. 
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INTRODUCTION 

The  purpose  of  this  note  is  to  outline  a  recently  developed  method 
for  formulating  and  solving  equations  of  motion  of  multibody  dynamic 
systems  that  is  well  suited  for  high  speed  dynamic  simulation  using 
parallel  processors.  An  outline  of  the  approach  is  given  here,  with 
references  cited  to  papers  that  develop  the  mathematical  foundations. 
Examples  involving  a  rigid  body  vehicle  system  real-time  simulation,  a 
geometrically  nonlinear  rotating  blade,  and  a  space  manipulator  are  used 
to  illustrate  application  of  the  method  and  to  indicate  computational 
efficiency  that  can  be  gained. 

OUTLINE  OF  METHOD 

References  1  and  2  present  a  recursive  algorithm  for  dvnamics  of 
rigid  body  systems,  either  open  loop  or  closed  loop,  based  on  the 
variational  form  of  multibody  equations  of  motion.  In  most  recent 
formulations  of  equations  of  motion  in  terms  of  joint  relative 
coordinates  f3],  equations  of  motion  of  bodies  in  a  mechanical  system 
are  first  written  in  terms  of  absolute  coordinates.  Kinematic 
relationships  between  absolute  velocity,  acceleration,  and  virtual 
displacement  and  corresponding  variables  in  terms  of  relative 
coordinates  are  applied  to  the  graph  theoretic  model  of  the  system,  to 
reduce  absolute  coordinate  equations  of  motion  to  coupled  equations  of 
motion  in  all  relative  coordinates  and  base  body  absolute  coordinates. 

In  Ref.  1,  Bae  observed  that  as  this  reduction  proceeds  from  tree 
end  bodies  toward  the  base  body,  coefficients  of  relative  coordinate 
variations  must  be  zero,  leading  to  equations  that  can  be  used  to  solve 
for  associated  relative  coordinate  accelerations  as  functions  of  inboard 
body  absolute  accelerations  and  velocities.  This  process  is  repeated, 
recursively,  from  tree  end  bodies  to  the  base  body,  to  obtain  a  reduced 
set  of  equations  of  motion  associated  with  the  base  body.  Upon  solving 
these  equations  of  motion  for  the  base  body  accelerations,  relative 
coordinate  accelerations  are  recovered,  moving  from  the  base  body 
outward  toward  tree  end  bodies.  Following  integration  of  these  relative 
coordinate  accelerations,  position  and  velocity  analysis  is  carried  out 


by  moving  from  the  base  body  outward  toward  tree  end  bodies,  completing 
one  step  of  the  system  equation  formulation  and  integration  process. 
Calculations  associated  with  this  recursive  process  are  implemented  with 
the  aid  of  symbolic  processing  software  to  yield  an  order  of  magnitude 
speed  up  over  absolute  coordinate  methods  [1,2].  Equally  important,  the 
tree  structure  associated  with  this  recursive  algorithm  leads  to 
definition  of  parallel  computation  tasks  that  can  be  carried  out 
effectively  on  a  parallel  processor  f 4 ] .  This  algorithm  is  currently 
yielding  to  real-time  vehicle  dynamic  simulation. 


A  RIGID  BODY  VEHICLE  EXAMPLE 

As  an  illustration  of  the  foregoing  approach,  consider  the 
schematic  representation  of  a  four  wheel  vehicle  with  double-A  arm 
suspension  subsystems  at  all  four  corners  and  rack  and  nin.ion  steering 
mechanisms,  as  shown  in  Fig.  1.  ^his  model  represents  the  Armv's  new 
HMMWV  multipurpose  vehicle.  It  consists  of  i  total  of  IS  Bodies,  22 
kinematic  joints,  anti  'losed  kinematic  loops.  Following  i  process  of 
cutting  selected  joints  to  open  kinematic  loops  [2],  a  tree  structure 
model  of  the  vehicle  is  obtained,  is  shown  in  Ei"'.  2.  This  model  leads 
to  10  parallel  paths  for  concurrent  processor  implementation. 


Figure  1  Schematic  Penresentat ion  Figure  2  Tree  Stucture 
of  a  Vehicle  Corresponding  to  Vehicle 

To  illustrate  parallel  processor  allocation,  the  parallel  task 
graph  of  Fig.  3  is  provided  [ 4 ] .  This  model  and  algorithm  is  currently 
running  on  an  Alliant  FX8  mini  supercomputer,  including  all  nonlinear 
effects,  with  approximately  12  millisecond  timesteps.  This  temporal 
resolution  is  more  than  adequate  for  real-time  simulation  of  vehicle 
systems,  accounting  for  approximate! v  50  Hertz  frequenev  response. 


A  GEOMETRIC  NONLINEAR  FLEXIBLE  BODY  EXAMPLE 

The  algorithm  outlined  above  for  rigid  body  applications  has 
recently  been  extended  to  multif lexible  body  systems  [5].  Deformation 
modal  coordinates  are  employed  in  this  formulation  to  account  for  linear 
elastic  deformation  effects,  relative  to  a  moving  reference  frame. 
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Figure  3  Parallel  Task  Graph  of  Recursive  Algorithm  for  Vehicle 

Nonlinear  coupling  due  Co  elastic  deformation  and  gross  motion  of  a  body 
reference  frame  are  accounted  for  in  the  formulation.  The  second  basic 
approach  outlined  above  for  rigid  bodv  systems  is  applicable,  with  the 
exception  that  deformation  mode  coordinates  are  eliminated  in  the  same 
manner  as  relative  coordinates  in  the  case  of  rigid  body  dynamics.  To 
illustrate  use  of  the  method,  a  long  slender  bar  shown  in  Fig.  A  is 
constrained  to  rotate  above  the  global  z  axis.  Its  motion  is  begun  at 
rest  in  the  configuration  shown  in  Fig.  4.  The  shaft  that  rotates  about 
the  z  axis  is  accelerated  to  a  steady  state  angular  velocity  of  four 
radiants  per  second  in  a  period  of  15  seconds.  Numerical  results  for 
tip  deflection,  relative  to  the  undeformed  state  of  the  bar,  are  shown 
in  Fig.  5  for  several  flexible  body  models. 


Tf  the  bar  is  modeled  as  a  single  body,  results  labeled  "linear 
model"  shown  in  Fig.  5  indicate  divergence.  As  two,  four,  and  six 
substructures  are  used  in  modeling  t^e  spinning  bar,  results  that 
account  for  geometric  stiffening  are  shown  in  Fig.  5.  In  these  later 
models,  each  substructure  is  assigned  a  bodv  reference  frame  that  moves 
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with  the  substructure,  approximating  an  Eulerian  reference  frame  that 
moves  with  the  system,  hence  permitting  capture  of  geometric  nonlinear 
effects.  These  results  agree  well  with  classical  theorv  and  numerical 
results  obtained  with  the  absolute  coordinate  PAPS  computer  code.  It  is 
interesting  to  note  that  even  though  there  are  no  parallel  tasks  to  be 
exploited  in  this  application,  CPU  time  required  for  the  recursive 
algorithm  in  the  six  substructure  model  is  one-thirteenth  that  required 
for  an  absolute  coordinate  approach.  This  suggests  that  not  on1,  v  can 
high  speed  simulation  of  flexible  body  systems  be  achieved,  but 
geometric  nonlinear  effects  can  be  captured  in  the  computational 
approach . 

A  FLEXIBLE  MULTIBODY  MANIPULATOR  EXAMPLE 

As  a  final  example  of  interest  in  space  structure  applications, 
consider  the  multi  flexible  body  model  of  a  shuttle  manipulator  shown  in 
Fig.  b  [4],  In  this  application,  the  two  long  intermediate  hcdies  are 
flexible.  The  rigid  base  rotates  about  an  axis  fixed  on  the  scuttle  and 
the  end  effector  is  modeled  as  a  rigid  body.  There  are  thus  four  movin' 
bodies  in  the  svsten,  each  connected  to  its  predecessor  hodv  bv  a 
revolute  joint.  As  a  simple  computational  illustration  [5]  relative 
position  actual  ^  r  s  at  the  Joints  wht»r  a  sped*' led.  and  dovl  ad  d  the 
flexible  body  simulation  compared  with  those  of  a  rigid  body  simulation 
in  Fig.  I.  As  noted,  substantial  elastic  deformation  effects  occur,  due 
to  the  gross  motion  and  nonlinear  coupling  effects  in  tKo  "’odel  . 


Figure  6  A  Flexible  Manipulator  Figure  7  Tip  Deviation  of 
with  Angle  Drivers  Flexible  Manipulator  System 

The  real-time  simulation  capability  currently  functioning  for 
realistic  vehicle  models,  coupled  with  successful  initial  demonstration 
of  capability  to  account  for  multi  flexible  body  effects,  suggest  that 
high  speed  and  perhaps  real-time  mu  1 1 i f 1 exible  body  simulation  will  he 
possible  in  the  near  future,  through  exploitation  of  emerging  parallel 
processor  architectures. 
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INTRODUCTION 

The  equations  for  the  motion  of  a  flexible  spacecraft  in  space  can 
be  described  conveniently  in  terms  of  quasi -coordinates .  The  nonlinear 
equations  can  be  treated  by  a  perturbation  approach,  resulting  : n  a  low- 
order  nonlinear  problem  for  the  "rigid-body”  maneuvering  and  a  .'.near 
time-varying  problem  for  angular  perturbations  from  the  rigid-body 
maneuver  and  elastic  vibration. 

THE  EQUATIONS  OF  MOTION  AND  THE  PERTURBATION  APPROACH 

The  motion  of  a  flexible  spacecraft  can  be  described  in  terms  of 
the  rigid-body  translation  and  rotation  relative  to  the  inertial  space 
of  a  set  of  body  axes  xyz  embedded  in  the  underformed  body  and  the 
elastic  motion  relative  to  the  body  axes.  Assuming  that  the  rigid-body 
translation  can  be  ignored,  the  equations  of  motion  can  be  written  in 
the  form  of  a  set  of  hybrid  (ordinary  and  partial  differential) 
Lagrange's  equations  in  terms  of  quasi-coordinates  [1] 

uJw  -  u)J"^o(r  ♦  u)v  dD  +  M  (la) 

(lb) 

where  oo  =  angular  velocity  vector,  v  =  elastic  velocity  vector,  r  = 

position  vector  of  a  point  in  the  body,  u  =  elastic  displacement  vector, 

p  =  mass  density,  M  =  torque  vector,  U  =  force  density  vector,  L  -  d i f — 

c  ~  -  T  .  ~~T 

ferential  operator  matrix  and  J  =  JDo(r  +  u)(r  +  u)  dD,  Jv  =  <Rp(rv  + 

vr^)dD.  Note  that  symbols  with  an  overtilde  represent  skew  symmetric 
matrix  counterparts  of  the  corresponding  vectors.  Moreover,  all 
quantities  in  Eqs.  (1)  are  in  terms  of  components  about  the  body  axes. 

Equations  (1)  can  be  conveniently  used  to  describe  the  maneuvering 
of  flexible  spacecraft.  The  equations  are  not  only  hybrid  but  also  non¬ 
linear.  It  should  be  noted,  however,  that  an  ideal  maneuver  is  one  in 
which  the  spacecraft  is  reoriented  as  if  it  were  rigid.  Hence,  a 
reasonable  assumption  is  that  the  rigid-body  rotations  are  generally 
large  and  the  elastic  displacements  are  generally  small.  This 


Jw  +  J_o(r  +  u)v  dD  =  -J  w  -  f^pvv  dD  - 
~JD  -  v~JD~ 

-  -  X •  •  -T 

p(r  +  u)  u)  +  ov  =  -ovoj  -  Cu  +  U 
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assumption  permits  a  perturbation  approach  [2)  to  the  solution,  in  which 
terms  involving  rigid-body  motions  are  regarded  as  being  of  zero-order 
in  magnitude  and  terms  involving  elastic  motions  are  regarded  as  being 
of  first-order  in  magnitude.  Hence,  we  assume  that 


“  =  “o  +  ~1’  J  =  J0  +  Jl’  -  =  ~0  +  ~1  (2a,b,c) 

where  J-  =  ("_ prr^dD,  J,  =  J"  p(ru^"  +  ur^)dD  =  J  ;  the  subscripts  indicate 
0  D  I  D  «  u 

the  order  of  magnitude.  Moreover,  u,  v,  U  and  Jv  are  small.  Intro¬ 
ducing  Eqs.  (2)  into  Eqs.  (1),  we  obtain  the  zero-order  equations 


J0~0 


=  ’“oVo  * 


*0 
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which  are  recognized  as  Euler's  equations,  and  the  first-order  equations 
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Equations  (3)  and  (A)  permit  design  of  a  maneuver  strategy  describ¬ 
ed  by  a  low-order  set  of  nonlinear  ordinary  differential  equations  and  a 
linear  set  of  hybrid  equations.  According  to  this,  Eq.  (3)  represents 
the  rigid-body  maneuvering  equations.  The  rigid-body  maneuvering  of 
spacecraft  has  been  discussed  extensively  in  the  literature  (see,  for 
example  [3]).  On  the  other  hand,  Eqs.  (A)  govern  angular  perturbations 
from  the  rigid-body  maneuver  and  elastic  displacements.  They  represent 
a  regulator  problem. 


As  an  illustration,  we  consider  the  single-axis  maneuver  of  the 
spacecraft  shown  in  Fig.  1.  The  pertinent  quantities  are  as  follows: 


So 
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Inserting  Eqs.  (5)  into  Eq.  (3)  and  introducing  an  obvious  kinematical 
relation,  we  obtain  the  maneuver  equations 


JxxC(t)  =  M  (t),  0  (t)  =  G(t)  (6a, b) 

where  9^(0  is  the  rigid-body  maneuver  angle.  Moreover,  inserting  Eqs. 
(3)  into  Eqs.  (A)  and  adding  appropriate  kinematical  relations,  we 
obtain  the  regulator  equations 
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pXli)  *  ov  =  /  u  -  PXft  +  U  ,  —  D  XL)  +  0V  =  4.  u  +  11 
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where  9,(t)  =  (9  (t)  9  (t)  9  (t)l  is  a  vector  of  perturbation  angles. 

“*1  x  y  z 

We  observe  that  in  this  particular  case  Eqs.  (7a),  (7b)  and  (7f)  are 
independent  of  Eqs.  (7c),  ( 7  d )  and  (7e).  This  impLies  that  the  varia¬ 
bles  u)  and  u  can  be  controlled  independently  of  the  variables  u>  ,  w 
z  y  x  y 

and  uz,  which  simplifies  the  control  problem  significantly. 

Equations  (6)  permit  a  maneuver  strategy  as  if  the  spacecraft  were 
rigid.  For  a  minimum-time  maneuver,  the  control  is  bang-bang  [4).  As  a 
result  of  selecting  a  rigid-body  maneuver,  ft  and  ft  can  be  regarded  as 
known  functions  of  time.  As  can  be  observed  from  Eqs.  (7),  the  effect 
of  ft  and  ft  is  to  introduce  time-dependent  coefficients  into  the  regula¬ 
tor  equations.  Hence,  the  net  effect  of  the  perturbation  approach  is  to 
replace  a  hybrid  nonlinear  problem  by  an  ordinary  second-order  nonlinear 
problem  and  a  hybrid  time-varying  linear  problem.  It  is  clear  that  the 
latter  alternative  is  much  more  advantageous. 

It  should  be  pointed  out  that  in  both  cases  described  above  spatial 
discretization  is  necessary  to  transform  a  hybrid  system  into  a  system 
of  ordinary  differential  equations.  Then,  the  question  arises  as  to  the 
nature  of  approximation  and  in  particular  whether  controls  designed  for 
a  discretized  system  are  capable  of  controlling  the  actual  structure. 
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The  cor:  epr  of  adaptive  structures  is  brought  up  in  connect  ion  wi 
the  r.eed  in.  tra  i  i  ght  we  i  ght  structural  systems  which  maintain  do.s;:v 
proper  t  ie  >  i : .  1  cor.  t  i, vsr.it  ions  without  human  intervention  when  subject*- 
to  dynamic,  thermal,  and  other  envi  rorur.ental  forces.  Examples  are  lar 
mtenna  s : rue t  tires  and  flexible  robotic  structures.  In  the  both  cases 
such  adapt  tvitv  would  allow  less  massive  structural  members  to  be 
employed  under  normal  loading  conditions  During  special  circumstance 
when  unusually  large  loads  an  encountered,  temporary  stiffening,  would 
allow  the  use  ot  less  st  irdv  structures,  resulting  in  large  savings  in 
their  cost,  and  in  increasing  their  mobility  and  efficiency. 

Within  the  framework  of  a  t inite -dimension.}  1  represent  at  ion  of 
structural  dynamics  the  adapt ivi tv  can  he  implemented  bv  the  dependent- 
of  -  he  stiffness  matrix  :  k  I  upon  f  he  ,-sp,-.  t,-d  load  .O',  or  expected 
<  programmed  ;  changes  in  conf i gurat ions .  i  e 

I  k  1  -  (  k  t  t )  ) 

where  the  dependence  upon  time  is  programmed,  in  advance.  In  order  to 
sustain  unexpected  loads  the  adaptive  structure  can  he  provided  bv 
feedback  force  control,  or  bv  a  parametrical  stiffness  control.  In  th 
last  case  instead  of  t 1 > 

1  k  1  ■  :  k  1 1 .  X  p  1  ’.  I 

in  which  X ^  is  a  generalized  coordinate 

A  variable  stiffness  ot  a  structural  member  can  be  achieved  bv 
appropriate  changes  in  pre-stress.  Indeed,  as  follows  from  finite 
e last ic i t v : 
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in  which  E  and  G  are  Che  Young's  and  the  shear  modulii,  respectivity 
and  Txx  is  the  normal  pre-stress,  etc. 

Hence,  the  stiffness  of  an  elastic  member  is  increased  by 
increasing  the  pre-stress  Txx.  This  can  be  implemented  by  cables, 
hvdraulic  actuators,  piezoelectric  actuators,  or  induced  thermal 
gradients.  A  variable  shape  of  a  structural  member  (for  instance,  a 
snake- like  robotic  arm)  is  achieved  by  an  appropriate  combination  of  a 
pre- tension  and  bending  moments  which  can  be  implemented  bv 
piezoelectrical  actuators  (for  fast  changes)  and  bv  thermal  gradients 
i  for  slow  changes). 

Since  adaptive  structures  are  characterized  by  t ime - dependent 
stiffness,  new  phvsical  phenomena  (such  as  a  parametrical  resonance i 
should  be  expected.  In  addition,  because  ot  time- dependent 
coefficients  in  the  governing  equations  for  transient  dynamics  all  the 
methodology  associated  with  the  modal  approach  as  well  as  the  stability 
ir.nl  vs  is  should  be  revised. 

In  application  to  robotic  arms  the  advantages  of  structural 
adaptivity  become  obvious.  Indeed,  a  rigid  segment  arm  can  operate 
with  loads  which  are  several  times  lighter  than  the  arm  itself,  while  a 
human  arm  (which  is  an  example  of  an  adaptive  "structure")  operates 
with  loads  several  times  heavier  than  its  own  weight.  One  of  the 
specific  problems  which  arises  in  connection  with  a  flexible  arm  is  the 
interaction  between  the  rigid  motion  and  elastic  micromotions.  Such  an 
interaction  leads,  for  instance,  to  a  false  impulse  of  the  arm  tip  due 
to  i'nite  time  of  the  elastic  wave  propagation  of  the  torque  applied  at 
the  arm  root.  (This  undesirable  effect  can  be  minimized  bv  introducing 
a  pre-stiffening  of  the  arm  before  rhe  maneuver.)  From  the 
mathematical  viewpoint  rigid  motions  are  governed  bv  the  Lagrange 
equations,  while  the  elastic  motion  described  by  appropriate  continuous 
models  (like  the  Timoshenko  beam,  or  the  Kirchhoff's  thin  rod)  are 
governed  by  partial  differential  equations.  (Finite  dimensional 
approximations  for  the  elastic  motions  as  well  as  the  Bernoul 1 i - Euler 
model  may  lead  to  a  loss  of  the  wave  propagation  effect.)  Due  to 
interactions  between  rigid  and  elastic  motions  the  Lagrange  ordinary 
differential  equations  and  partial  differential  equations  are  coupled: 
the  Lagrange  equations  incude  an  additional  (reactive)  force  depending 
upon  the  integral  effect  of  the  elastic  motions,  while  the  partial 
differential  equations  parametrically  depend  on  the  angular  velocities 
of  the  rigid  motions.  Since  the  displacements  due  to  the  rigid  motions 
are  not  small,  the  total  system  is  non-linear. 

Traditionally  the  robotic  research  is  being  performed  bv  the 
community  of  electrical  engineers.  However,  even  a  brief  review  of 
flexible  robotic  structures  demonstrates  the  existence  of  complex 
problems  which  require  deep  understanding  of  structural  dynamics  and 
careful  selection  of  appropriate  mathematical  models.  That  is  whv  now 
is  probably  a  good  time  for  the  community  of  structural  engineers  to 
offer  their  assistance  in  the  development  of  a  new  generation  of  robots 
with  adaptive  flexible  arms. 
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Two  types  of  specific  electromechanical  actuators  are  described  and 
discussed.  A  proof  mass  actuator  and  an  electric  motor  are  examined 
in  terms  of  the  amount  of  damping  each  produces  in  a  specific 
structural  control  experiment.  Theoretical  and  experimental  values  of  actuator 
produced  damping  are  examined.  The  effects  of  actuator  dynamics  on  control 
law  implementation  are  noted.  In  addition,  a  theoretical  parameter 
study  of  the  dynamic  response  of  a  DC  motor  controlling  a  flexible 
model  of  a  beam  are  summarized. 

The  proof  mass  actuator  hardware  consists  of  a  microprocessor  controlled, 
solenoid  like  mass  with  displacement  and  acceleration  sensors. 
Local  digital  control  is  used  to  cause  the  actuator  to  provide  velocity  feedback 
at  a  point.  This  form  of  control  is  capable  of  providing  roughtly  a  10%  increase 
in  damping  for  various  nodes.  Transfer  function  analysis  of  the  actuator  revails 
dynamics  of  the  form 

F(s)  G1G2mrs2 


V(s)  mrs2+G1G2C-|S  +  G1G2k 

where  F  is  the  force  output,  V  is  the  input  voltage  and  G-| ,  G2  Mr.  c1  and  k  are  the 

various  physical  parameters  of  the  actuator  and  amplifier  [1],  The  actuator 
dynamics  become  important  if  the  first  natural  frequency  of  the  structure  falls  below 
the  break  frequency  of  the  actuator  (about  2.5  hertz). 

A  similar  analysis  is  also  presented  for  an  electric  motor  used  as  a  hinge  and 
controlled  for  vibration  suppression  during  slewing  maneuvers.  The  actuator 
consists  of  a  DC  permanent  magnet  motor  and  tachometer  of  approximately  16% 


of  the  mass  of  the  structure  [2],  The  electric  hinge  is  located  at  the  8  foot 
mark  of  a  12  foot  long  flexible  aluminum  beam.  The  motor  and  gear 
combination  is  first  compensated  for  frictional  losses,  then  active  structural 
control  was  accomplished  by  viscous  damping  introduced  through  simple 
tachometer  feedback.  Damping  increases  as  much  as  a  factor  of  3  (from 
8.3%  to  23%)  were  measured  [3]. 

A  very  simple  model  of  a  motor  controlling  a  load  is  used  to  illustrate  how 
damping  is  introduced  by  the  motor  into  the  load  dynamics.  The  damping 
introduced  in  this  fashion  constitutes  a  structure,  control  law  and  actuator 
interaction,  as  the  equivalent  viscous  damping  parameter  depends  on 
coefficients  characterizing  the  load  a  well  as  the  motor.  A  more  detailed  model 
of  the  motor  illustrates  that  when  the  load  inertia  is  comparable  to,  or  greater 
than,  the  motor  inertia,  substantial  actuator  structure  interaction  occurs  greatly 
effecting  control  law  performance.  The  details  of  this  interaction  depend  on  all 
of  the  systems  parameters  and  hence  are  not  simple.  State  feedback  optimal 
control  is  used  to  determine  structure  and  control  parameter  that  yields 
minimum  power  usage  for  constant  speed  of  response.  The  procedure  has 
been  repeated  for  a  flexible  load  [4], 
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Comments  in  response  to  a  questions  raised  during  the  presentation  (Crawley) 

A  q  .estion  was  raised  regarding  the  nature  of  the  high  gain  instability  of  a  proof  mass 
actuator. 

As  was  noted  in  the  lecture,  the  proof  mass  actuator  illustrates  a  high  gam  instability. 
That  is,  as  the  gain  of  the  actuator  is  increased,  the  stability  of  the  structure  actuator 
system  is  lost.  This  phenomena  is  explained  in  terms  of  the  interaction  between  the 
control  and  the  structure.  The  actuator  can  be  viewed  simplistically  as  an  added 
degree  of  freedom  which  "absorbs"  vibration  of  the  structure  modes  by  adding 
damping  to  these  modes.  As  the  control  gain  is  increased,  more  of  the  damping  force 
of  the  actuator  is  transferred  to  the  structural  modes.  This  happens  at  the  reduction  of 
the  damping  force  of  the  actuator  mode.  Eventually  the  actuator  mode  has  no  damping 
force  and  the  system  (structure  plus  actuator)  beomes  unstable  (actually,  marginally 
stable).  If  the  gain  is  increased  further  the  actuator  modal  damping  becomes  negative, 
causing  an  instability  in  the  system  [1],  This  is  predicted  analytically  and  observed 
experimentally. 

(1)]  Zimmerman,  D.C.  and  Inman,  D.  J.,  "Actuator  Structure  Ineraction.  "Proceedings 
of  the  Second  NASA/DOD  Control/Structure  Interaction  Technology  Conference, 
November  1987. 
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INTRODUCTION 

Most  active  controllers  developed  to  control  large  structures  are 
subject  to  constraints.  For  instance,  control  characteristics  of 
proof-mass  actuators  are  dominated  by  the  nature  of  the  constraints. 

To  find  the  optimal  or  suboptimal  control  laws  for  such  controllers  can 
be  a  formidable  task.  Intuition  of  a  designer  plays  an  importruit  role 
[1]  and  the  design  method  may  vary  dramatically  in  accordance  with  the 
constraints  of  the  system.  Therefore,  one  may  wish  to  have  a 
systematic  methodology  to  solve  the  control  problems  subject  to  control 
force  and  state  variable  constraints.  Limiting-performance  / 
minimum-time  (LP/MT)  control  calculates  the  optimal  control  force  as  a 
function  of  time  for  a  known  system  with  initial  conditions,  subject  to 
certain  constraints  and  external  excitations  while  minimizing  a  given 
performance  index.  Although  the  LP/MT  control  gives  optimal  open  loop 
control  for  the  system,  it  is  desirable  to  develop  a  closed  loop 
control  which  has  more  practical  value.  Perhaps  a  system 
identification  technique  can  be  used  to  establish  a  suboptimal  feedback 
control  law  based  on  the  LP/MT  control  characteristics. 


LIMITINC-PERFORMANCE/MINIMUM-TIME  CONTROL 

A  brief  description  of  the  problem  which  the  method  can  solve  is 
given  here.  Further  information  on  the  method  can  be  found  in  [2]  and 
[3].  A  linear  vibrating  system  with  n  degrees  of  freedom  subject  to 
arbitrary  external  excitations  f(t)  and  control  forces  u(t)  is 
expressed  in  the  first  order  system  of  differential  equations 

s(t)  =  As( t)  +  Bu( t )  +  Cf(t)  (1) 

where  s(t)  is  an  n-dimensional  state  vector.  A,  B,  and  C  are  n  x  n,  n  x 
nu  and  n  x  nf  constant  coefficient  matrices.  The  quantities  nu  and  nf 
are  the  number  of  control  forces  and  excitations,  respectively. 

Constraints  are  imposed  on  the  dynamic  system  under  study.  The 
format  of  the  constraints  is 

+  02—  +  Q3I  £  Y[j  for  tQ  <  t  <  tj.  (2) 

where  y^  and  y^  are  nc-dimensional  lower  and  upper  constraint  vectors. 
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Qj ,  Q9  and  are  nc  x  n.  nc  x  nu  and  nc  x  nf  constant  coefficient 

matrices,  ;ind  c  and  t  .  are  given  initial  and  final  time, 
o  t 

The  LP/MT  control  finds  an  optimal  control  u( t )  which  will 

transfer  an  initial  state  s(t  )  =  s  to  a  desired  final  state  s(t  )  = 

-v  o  o  -  f 

s j.  in  the  minimum  time  while  extremizing  :i  given  performance  index  of 
the  form 

max  .y.  y  y 

Minimize  J  =  {t  <  t  <  t  j.  Ig^s  +  +■  p^f  | }  (3) 

where  B3  are  given  n,  nu  :ind  nf  constant  coefficient 

vectors . 


SYSTEM  IDENTIFICATION  OF  SUBOFI I  MAI .  FEEDBACK  CONTROL  PARAMETERS 

Since  the  l.P/MT  control  gives  the  best  possible  or  "limiting" 
response  of  a  system,  it  would  appear  to  be  reasonable  to  base  a 
control  system  on  the  LP/MT  control  characteristics.  However,  due  to 
uncertainties  in  control  problems,  open  loop  control  such  as  the  LP/MT 
control  may  not  be  applicable  in  practice,  unless  real-time  computing 
power  for  the  LP/MT  control  is  available.  To  overcome  this  difficulty, 
parttmeter  identification  to  find  suboptimal  feedback  control  laws  based 
on  the  LP/MT  control  characteristics  is  possible.  Note  that  the 
prescribed  constraints  have  been  taken  into  account  when  the  optimal 
LP/MT  control  is  computed. 

Consider  a  linear,  time  invariant  dynamic  system  represented  by  a 
set  of  difference  equations 

s(k  +  l)  =  Gs(k)  +  Hu(k)  (T) 


where  s(k)  and  u(k)  are  vectors  of  dimensions  n  and  m,  respectively;  G 
and  H  are  constant  matrices  of  appropriate  dimensions. 

From  the  LP/MT  control  characteristics,  the  optimal  time  responses 
s  (k)  and  u  (k)  are  obtained.  Consider  a  linear  controller, 


u(k)  =  Ks(k) 


(5) 


where  K  is  an  m  x  n  feedback  gain  matrix. 

Since  m  controllers  are  considered  and  optimal  control  forces  are 
available  for  each  controller,  it  is  possible  to  proceed  controller  by 

controller.  For  controller  i  (i  =  1,2 . m),  a  suboptimal  linear 

control  law  to  be  determined  is  described  by 
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for  controller  i 

Us*(k) 

=  k. 
-] 

s*(k) 

(7) 

k.  =  i1^  row  of  K  matrix 

- 1 


Here  it  is  supposed  that  the  controller  i  is  constrained  as 
(See  Eq.  (2)) 

Umin.  <  u.(k)  <  Umax .  (S) 

l  -  1  '  -  i 

Define 

re.(k)  =  u*(k)  -  Us.(k)  (9) 

* 

where  u.(k)  is  an  optimal  control  force  for  controller  i. 

Then  a  constant  feedback  gain  matrix  K  in  Eq.(5)  is  selected  to 
minimize 

N-l  9 

RE.  =2  [re. (k)],  i  =  1,2 . m  (19) 

1  k=l  1 

This  identification  method  handles  problems  efficiently,  since 
once  the  optimal  LP/MT  control  characteristics  are  known,  it  is  not 
necessary  to  perform  a  complete  structural  analysis.  It  is  expected 
that  this  parameter  identification  technique  will  give  a  stable 
feedback  control  law  for  a  large  structural  system.  To  ensure 
stability  the  gain  matrix  K  will  be  constrained  when  a  least  squares 
fit  is  performed.  Also,  to  increase  the  fitness  nonlinear  gain 
matrices  may  be  superimposed  on  a  linear  gain  matrix. 
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Notes  on  Presentations  and  Discussions  in  the 
Modelling  Session 

Contributed  by 
Edward  J .  Haug 

P.M.  BAINUM 

Peter  Bainun  summarized  the  results  from  two  preceding  papers,  as 
indicated  in  the  abstract  contained  in  the  proceedings. 

During  discussion,  Meirovitch  inquired  as  to  what  rigid  body  nodes  are 
accounted  for.  Bainun  implied  only  rotational  rigid  body  nodes  were  included, 
since  only  they  have  significant  effect. 

Bryson  questioned  the  speaker  regarding  interpretation  of  poles  and  the 
significance  of  high  frequency  nodes  in  the  analysis.  Bainun  responded  that 
the  scope  '.if  the  study  was  set  by  NASA-Langlev ,  the  sponsor. 

Eubanks  inquired  about  nonselfad jointness  of  the  model.  Bainum  responded 
that  he  was  not  looking  at  nonsel f ad i ointness .  Junkins  and  Meirovitch 
clarified  the  issue  regarding  self adiointness  in  open  loop  control. 

Bryson  extended  his  earlier  comments  regarding  structural  frequency 
ranges  of  interest  and  their  coupling  with  linear  effects.  He  argued  that  the 
degree  of  coupling  at  higher  freouencv  ranges  is  insignificant.  Bainum  agreed 
that  high  frequency  effects  are  essentially  linear  in  structural  response. 

Amos  Inquired  about  the  significance  of  coupling  between  deformation  and 
rigid  body  motion.  Bainum  responded  that  the  application  being  considered  did 
not  concern  slewing  effects;  i.e.,  only  small  gross  motion  was  of  interest. 

J.H.  WILLIAMS,  JR. 

Jim  Williams  summarized  considerations  regarding  wave  propagation 
modeling  and  feedback  control,  as  compared  to  nodal  techniques.  He 
specifically  related  modeling  to  wave  length  and  structural  dimension 
parameters . 
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Veirovitch  challenged  Williams  on  how  to  excite  nodes  in  space.  Wi 1  liars 
responded  that  actuators  provide  a  significant  excitation  of  waves  and  nodes. 

Eubanks  challenged  Williams  regarding  classification  of  problems  vrith 
only  wave  speed  (wave  length.)  and  characteristic  structural  length.  He  argued 
that  parameters  associated  with  the  nature  of  load  application  were  equally 
important.  Crowley  reinforced  his  opinion,  observing  that  the  manner  of  load 
application  dictates  structural  response  and  the  approach  taken  to  modeling. 

Highland  argued  that  the  wave  approach  is  important  and  has  substantial 
potential  application  in  control  and  should  be  given  serious  consideration. 

He  voiced  the  opinion  that  the  wave  approach  is  not  being  given  fair 
cons iderat ion ,  in  comparison  with  wave  approaches. 

Meirovitch.  argued  that  dispersion  of  waves  is  important  and  does  not 
appear  to  be  fully  addressed.  Williams  responded  that  disnersion  is 
important,  but  nondisDersive  effects  are  also  of  considerable  importance. 

Haug  and  Craig  argued  that  deformation  shapes  other  than  global  vibration 
modes  need  to  be  given  consideration  for  modeling  structural  response.  Static 
correction  modes  used  to  characterize  local  deformation  fields  due  to  low 
transmission  across  structural  component  boundaries  are  significant  and  can 
lead  to  much  improved  modal  models,  as  compared  to  truly  normal  vibration 
modes . 


E.J.  HAL'G 

Ed  Haug  summarized  recent  developments  in  recursive  modeling  for  parallel 
processor  multi  flexible  body  application,  as  per  tKe  abstract  contained  in  tKe 
proceed i ngs . 

Von  Flotow  observed  that  Cain  and  co-workors  at  Stanford  have  observed 
geometric  stiffening  effects  outlined  in  the  presentation.  Haug  responded 
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Chat  the  purpose  of  presenting  these  results  was  two-fold.  First  to  counter 
an  earlier  assertion  that  multiflexible  body  codes  could  not  capture  geometric 
nonlinear  effects  and  second  to  demonstrate  that  use  of  piecewise  F.ulerain 
reference  frames  with  elastic  structural  models  relative  to  these  frames  can 
account  for  geometric  nonlinear  effects,  in  a  unified  and  computationally 
effective  manner. 

Siljak  observed  that  a  similar  automotive  simulation  application  has  been 
used  with  individual  suspension  suhsvstems  for  high  speed  computational 
algorithms.  Haug  resnonded  that  the  algorithm  presented  has  led  to  an  order 
of  magnitude  sneed  uo  in  computation  and  an  ahilitv  to  account  for  ^uch  higher 
system  frequency  response.  He  agreed  that  alternate  physical  views  of 
mechanical  systems  can  lead  to  concurrent  computational  paths. 

The  question  was  raised  as  to  how  two  orders  of  magnitude  of  computa¬ 
tional  speedup  in  the  vehicle  simulation  reported  was  achieved.  Haug 
responded  that  nearly  an  order  of  magnitude  was  gained  through  the  recursive 
algorithm,  apart  from  parallelism.  An  additional  factor  of  approximately 
three  was  gained  through  parallel  computation.  The  remaining  speedup  is 
associated  with  higher  speed  processors,  which  are  expected  to  provide  still 
further  gains  in  the  near  future. 

L.  ME IRQVITCH 

Meirovitch  expanded  upon  the  quasi-coordinate  formulation  for  maneuvering 
equations  for  spacecraft  that  is  presented  in  the  abstract  of  these 
proceedings . 

Dwyer  asked  if  he  could  start  with  a  finite  element  modal  characteriza¬ 
tion  of  a  large  scale  finite  element  model  of  a  structural  component  and 
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recover  coupling  effects.  Meirovitch  indicated  that  he  has  not  had  experience 
with  such  applications,  but  felt  that  it  was  possible. 

Venkavva  asked  if  the  proportion  of  pass  in  the  hub  of  the  spacecraft 
influences  the  equations  of  motion.  Meirovitch  indicated  that  a  more 
important  effect  is  the  regulator  that  couples  gross  notion  with  deformation 
effects . 

W.  PILKEY 

Pilkev  made  a  presentation  on  a  linear  programming  approach  to  open  loon 
controller  design,  as  in  the  abstract  in  these  proceedings.  He  emphasized 
that  the  open  control  design  techniaue  being  presented  is  independent  of  the 
number  of  degrees  of  freedom  of  the  structural  model.  He  pointed  out  that 
well  established  ’inear  programming  methods  for  solution  can  he  applied 
effect ivelv. 

Junkins  inquired  as  to  whether  the  approach  treats  bang-bang  control. 
Pilkey  indicated  that  his  approach  smootKs  hang-bang  control  effects. 

Eubanks  asked  if  the  formulation  :s  linear.  Pilkey  indicated  that  it  is, 
unless  nonlinear  kinematic  effects  are  accounted  for. 

Junkins  pointed  out  that  pole  placement  optimization  yields  nonlinear 
constraints.  Pilkey  agreed. 
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Notes  on  Presentations  and  Discussions  in  the 
Modelling  Session 


Contributed  by 
John  L.  Junkins 


P.  Bainum  -  Orbital  ['Meets 


Tunc  0,  R  weight  matrices  to  place  closed  loop  poles.  Used  fixed  stability  maruin  1  all 
eigenvalues.  Look  at  effects  ot  neglecting  orbit  attitude/  vibration  coupling  effects.  Which 
effects  dominate  depends  upon  how  large  deflections  are. 

Q.  What  are  the  bandwidths  ot  actuators,  sensors,  and  controllers'.’ 

A.  Don't  know,  off  hand. 

Q.  How  can  you  design  feedback  controller  without  worrying  about  bandwidth  and  roll  ot! 
characteristics? 

A.  We  didn't  pretend  this  was  an  optimal  controller.  We  simply  assumed  a  tor  example 
control  law. 


J,  W'illiams  -  Wave  Approach 

Suggests  wave  approach  is  appropriate  when  wavelength  (  )  of  excitation  is  much  less  than 
characteristic  length  (L)  ot  the  structure. 

F:.  liana  -  Multi-bodv  Dynamics 

Presented  a  new  formulation  tor  chain  type  topologies  which  facilitates  parallel  processing. 
Reported  two  order  ot  magnitudes  speedup  in  computing  motion  ot  a  system  with  11)  parallel 
paths,  some  of  the  speedup  due  to  taster  processor. 


L.  Mcirovitch  -  Maneuvers  Quasi-coordinates 

Presented  a  version  ot  Lagrange's  equations  using  quasi-coordinates  and  a  perturbation 
method.  I  he  zeroth  order  solution  is  the  rniid  body  part,  the  flex  body  effects  enter  as  a  lust 
order  correction. 

Q.  Can  the  zeroth  order  equations  be  modified  to  include  certain  dominant  flex  body 
coupling  effects? 

A.  Yes' 


W.  D.  Pijkey  -  System  Identification  of  Suboptimal  Feedback  Gains 


Find  discretized  open  loop  control  subject  to  inequality  constraints  on  control  t'orccv  positio 
of  proof  mass,  response  of  the  system.  Curve  tit  typical  response  history  to  litul  linear 
feedback  gains  "consistent"  with  the  optimal  open  loop  solution. 

Q-  Your  open  loop  trajectory  is  a  function  of  the  initial  conditions;  how  do  you  account  for 
this  in  trying  to  curve  fit  to  determine  gains’ 

A.  A  family  of  typical  initial  conditions  are  considered 


M.  Zak  -  Adaptive  St r 1 1 c t \ i r c s 

Presented  an  idea,  not  results.  Structures  must  bear  peak  loads  temporarily.  Consider 
temporary  increase  in  stillness  to  thereby  save  weight.  He  developed  equations  ot  motion 
very  similar  to  Meirovitch's.  Flexible  motions  are  perturbations  ot  rieul  motions  Heads  to 
nonlinear,  time  varying  equations. 

Q-  The  basic  idea  is  not  new,  is  it? 

A.  No,  but  it's  time  we  did  something  with  the  idea 


Dan  Inman  -  Actuators 


Discussed  Haviland.  UVn's  proof  mass  actuators.  Actuator  won't  work  below  2.5  Hz  and  has 
some  high  frequency  instability  which  has  not  been  carefully  identified  (by  Inman).  He 
invites  suggestions  for  experiments. 
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INTRODUCTION 

The  AFOSR- sponsored  ISAAC  program  at  the  Aerospace  Corporation  has 
been  pursuing  issues  which  arise  in  the  simultaneous  design  of  struc¬ 
tures  and  control  systems  for  larqe  space  structures  using  a  mathema¬ 
tical  programming  code  (see,  for  example  [6]).  There  are  many  practical 
advantages  to  such  integrated  design,  such  as  tuning  the  structure  to 
directly  improve  the  closed  loop  performance  measure.  This  talk  will 
focus  on  key  elements  in  our  work,  and  in  particular  those  elements 
which  distinguish  our  emphasis  from  that  of  other  workers  in  the  field. 

The  idea  is  to  simultaneously  optimize  some  cost  criterion  and  sa¬ 
tisfy  various  design  constraints.  Other  researchers  have  pursued  simi¬ 
lar  approaches  (see,  for  example,  [4]  and  references  therein).  It  is 
relatively  easy  to  list  the  desired  qualitative  characteri sties  of  good 
designs.  Development  of  an  implementation  involves  translation  of  these 
characteristics  into  quantitative  measures,  and  then  the  choice  of  which 
characteri stic(s)  to  implement  as  the  cost  criterion  and  which  to  imple¬ 
ment  as  constraints.  The  next  logical  step  is  to  become  concerned  with 
numerical  issues:  ideas  involving  the  solution  of  mathematical  program¬ 
ming  problems  can  easily  be  very  attractive  in  concept  but  unworkable 
computationally. 

In  our  ISAAC  work  we  have  attempted  to  address  both  conceptual  and 
computational  issues.  The  following  are  issues  we  have  considered  which 
may  be  classified  as  conceptual:  (1)  definition  of  system  robustness 
and  other  performance  criteria,  (2)  structure,  complexity  and  parametri- 
zation  of  controller,  and  (3)  structure,  complexity  and  parametri zat i on 
of  physical  spacecraft  structural  model.  Numerical  issues  we  have  con¬ 
sidered  include:  (4)  implementation  of  the  constraint  of  closed-loop 
stability,  (5)  algorithms  for  updating  structural  models,  (6)  simplified 
continuum  models  for  complex  structures,  and  (7)  evaluation  of  con¬ 
straint  and  criterion  functions  and  their  derivatives.  In  our  presenta¬ 
tion  we  shall  elaborate  on  these  issues  and  indicate  what  we  think  could 
be  accomplished  in  the  future. 

SUMMARY  OF  ISSUES 

Definition  of  Criteria:  The  translation  of  design  goals  into  criterion 
functions  is  central  to  the  success  of  optimization-based  design. 

Disturbance  attenuation,  insensitivity  to  parameter  variation,  band¬ 
width,  saturation  bounds,  tracking  performance  and  stability  are  typical 
examples  of  desirable  characteri sties  or  limitations  from  the  control 
design  point  of  view.  All  of  these  can  be  implemented  in  a  reasonably 
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direct  manner  using  the  frequency-dependent  maximum  singular  values  of 
transfer  function  matrices  [5,7],  There  is,  for  example,  no  need  to  use 
indirect  measures  of  disturbance  attenuation  or  robustness  such  as  damp¬ 
ing  ratios  as  in  [4] . 

From  the  structural  designers'  point  of  view  it  is  necessary  to  meet 
constraints  involving  permissible  ranges  of  parameter  variations  due  to 
requirements  on  physical  dimensions,  loading,  mass.  etc.  These  con¬ 
straints  are  generally  easy  to  implement  in  an  obvious  way. 

Controller  Structure:  Quadratic  regulator  theory  has  been  used  in  other- 
work  (e.g.  [4,8]),  and  this  has  dictated  the  compensator  structure.  Our 
approach  using  the  Q-parametri zation  [9]  of  stabilizing  compensators 
allows  a  greater  variety  of  compensator  structures.  It  does  still  leave 
open  the  question  of  how  to  best  limit  compensator  dynamic  order.  [3] 
argues  that  in  a  digital  implementation  of  such  an  approacn  there  is  no 
practical  need  to  limit  order. 

Spacecraft  Structural  Model :  We  have  worked  both  with  continuum  and 
finite  element  models.  The  mathematical  programming  approach  involves 
imposing  constraints  on  structural  parameters  based  upon  physical  signi¬ 
ficance.  The  presence  of  a  multitude  of  parameters  in  a  preliminary 
design  poses  an  obstacle  both  to  understanding  and  computation.  The  use 
of  continuum  models  [1]  can  provide  a  more  ready  conceptual  understand¬ 
ing  of  the  structure,  an  easier  handle  on  the  important  design  parame¬ 
ters,  and  means  of  more  readily  obtaining  a  feasible  design  point. 

There  are  also  cases  in  which  a  reasonable  continuum  model  cannot  be 
found . 

Closed  Loop  Stability  Constraint:  This  is  necessary.  The  quadratic 
regulator  theory  ensures  closed  loop  stability,  but  it  gives  a  fixed 
(and  high  dynamic  order)  compensator  structure.  Other  approaches  have 
involve  checking  stability  conditions  at  each  iteration  step:  for  exam¬ 
ple,  testing  for  positivity  of  certain  matrices  [8],  The  Q-parameteri - 
zation  provides  stability  automatically,  so  that  the  constraint  is  im¬ 
plicit  in  the  compensator  structure. 

Structural  Model  Updates:  In  the  process  of  varying  structural  parame¬ 
ters  when  using  frequency  domain  criteria,  it  becomes  necessary  to  re¬ 
compute  the  frequency  response  of  the  whole  system.  There  are  various 
possibilities  for  both  approximate  and  exact  frequency  response  updates, 
with  or  without  modal  analysis.  We  have  investigated  an  approach  using 
exact  frequency  response  updates,  which  is  computationally  efficient  for 
localized  structural  parameter  variations.  Conceptually  similar  tech¬ 
niques  appear  in  the  literature  for  exact  eigenvalue  updates  with  local¬ 
ized  modi f ications  [2] . 

Continuum  Mode  1 s :  Complex  structures  may  require  the  interconnect  ion  of 
several  such  models,  perhaps  with  finite  elements  models  as  well.  Effi 
cient  generation  of  frequency  responses  for  such  models  would  be  very 
useful . 

Evaluation  of  Criteria:  Numerical  evaluation  of  constraint  and  crite¬ 
rion  functions  can  involve  large  computations  at  each  step  of  the  itera- 


tion,  such  as  eigenvalue  problems,  singular  value  decompositions  or  mul¬ 
tiple  matrix  inversions.  One  generally  also  needs  to  differentiate  the 
results  of  these  computations  with  respect  to  parameters. 

CONCLUSIONS 

There  is  obviously  more  to  gain  by  simultaneously  tuning  structural 
and  contoller  designs.  Closed  loop  performance  criteria  which  directly 
measure  desired  characteristics  are  available.  The  big  issues  are  para- 
metrization  of  controller  and  structure,  and  computational  techniques 
for  evaluation  of  criteria,  constraints,  and  their  derivatives. 
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During  our  presentation  at  the  forum,  a  number  of  concerns 
pertaining  to  our  methodology  and  design  example  were  expressed  by 

several  forum  participants.  In  this  summary  we  briefly  address  those 

concerns . 

Compensator  C omp 1  ex i t y 

We  noted  that  Q ■ paramet rizat ion  design  generally  result  .  m  ,i 

compensator  which  is  full  (each  input  affects  each  output)  and  has  order 
higher  than  the  design  model  of  the  structure.  Based  upon  the  work  of 
Boyd,  et.  al.[l],  we  made  the  comment  that  if  the  resulting  compensator 
could  be  constructed  as  a  digital  finite  impulse  response  (FIR)  filter, 
then  it  would  not  be  necessary  to  reduce  high  compensator  order.  This  is 
because  VLSI  FIR  chips  are  already  being  marketed  which  contain  large 
numbers  of  delay  taps  [1].  The  concern  was  raised,  however,  that  such 
implementations  have  not  been  qualified  for  spaceflight,  and  hence 
complex  compensators  would  not  be  acceptable. 

We  agree  that,  when  acceptable,  compensator  simplification  is 

desirable.  However,  for  extremely  high  performance  mission  requirements 
compensator  complexity  will  be  unavoidable  due  to  increasing  complexity 
of  the  design  model  of  the  structure.  By  the  time  that  most  missions 
requiring  precision  structural  control  are  ready  to  fly,  VLSI  FIR 
hardware  will  most  likely  be  space  qualified  due  to  requirements  of 
advanced  communications  payloads. 

T OYSAT  Des ign  Example 

We  showed  a  diagram  of  the  TOYSAI  moael  wnich  is  tne  Dasis  of  our 
development.  This  consisted  of  two  2-dimensional  four-bay  truss 
appendages  attached  to  a  central  body.  At  each  truss  tip  a  vertical 
position  sensor  is  located,  and  a  rotation  sensor  is  located  at  the 
upper  center  body.  Here,  vertical  and  horizontal  are  defined  in  the 
TOYSAT  reference  frame  with  the  horizontal  direction  in  parallel  with 
the  undeflected  truss  center  line.  Two  force  and  one  torque  actuator  are 
respectively  collocated  with  the  sensors.  As  presented,  the  design  goal 
was  to  minimize  truss  tip  vertical  deflections  and  center  body  rotation 
responding  to  a  vertically-oriented  stat i on - keepi ng  thruster  located  at 
the  bottom  of  the  center  body. 


It  was  pointed  out  during  the  forum  that  no  center  body  rotation 
would  be  induced  by  the  thruster  disturbance  due  to  the  position  and 
orientation  of  the  thruster. 

We  appreciate  that  this  point  was  brought  out.  Earlier  in  our  work 
the  structural  dynamics  members  of  our  research  team  chose  the 
disturbance  source  to  be  a  bang-bang  slew  torque.  In  the  process  of 
coding  the  design  this  was  inadvertently  taken  to  be  a  station-keeping 
thruster.  We  have  since  corrected  the  disturbance  source  in  the  design 
example  so  that  all  three  outputs  will  respond  to  a  center  body  slew 
torque. 

We  might  add  that  even  in  the  example  presented,  the  design 
procedure  would  not  be  invalidated  because  the  truss  tips  will  respond 
to  the  thruster.  The  disturbance  to  output  transfer  matrix  will  be  3  by 
1  with  one  element  equal  to  zero.  We  have  verified  this  by  looking  at 
our  current  computational  results.  The  singular  value  of  this  matrix 
will  be  non-zero,  and  hence  the  objective  function  will  reflect  the 
truss  tip  responses. 

Multiple  Ext r  em a 

We  emphasized  the  use  of  constrained  optimization  to  select 
compensator  and  structural  design  parameters.  The  question  of  the 
possible  existence  of  multiple  extrema  was  raised. 

At  this  point  we  have  no  analytical  results  governing  the  number  of 
extrema  and  which  are  maxima  or  minima.  The  existence  of  multiple 
extrema  is  a  common  problem  in  mathematical  programming.  In  most 
practical  situations  one  must  run  the  optimization  for  a  variety  of 
initialization  points  and  pick  the  optimum  from  these,  with  no  assurance 
of  having  found  a  global  optimum.  Engineering  insight  will  improve  the 
likelihood  that  a  chosen  initial  point  will  converge  to  the  global 
optimum.  The  Q-  parametrization  design  method  [2],  by  virtue  of 
providing  some  ability  to  fix  input-output  structure,  may  facilitate 
such  insight.  If  the  form  of  the  compensator  is  chosen  to  be  digital  FIR 

[1],  then  most  objectives  and  constraints  will  be  convex  in  at  least  the 
controller  parameters,  although  still  not  in  the  structural  parameters. 
Convex  programming  problems  have  unique  minima  [3], 
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Modal  descriptions  of  dynamic  behavior  are  the  basis  of  the 
majority  of  proposed  structural  control  designs.  While  modal  analysis 
is  a  useful  tool.  limitations  to  its  applicability  exist.  Modal 
analysis  of  dense  spectrum  systems.  typical  of  proposed  large  space 
structures,  is  known  to  be  extremely  sensitive  to  small  parameter 
perturbations.  Some  analyses  indicate  that  hundreds  of  modes  of  such 
space  structures  will  contribute  to  performance  degradation.  Many  of 
these  modes  are  considerably  beyond  the  range  where  they  may  be 
confidently  modelled  using  this  technique.  Ignorance  of  the  eigen 
properties  of  these  modes  can  lead  to  unstable  control-structure 
interaction. 

Recent  work  performed  at  the  M.I.T.  Space  Systems  Laboratory1  has 
dealt  with  control  design  based  on  wave  propagation  models  of  flexible 
structures.  Reflection  and  transmission  properties  of  performance 
critical  locations  are  actively  altered  in  order  to  meet  mission 
requirements  regarding  dynamic  isolation,  energy  shunting,  and  vibration 
suppression.  Using  an  input/output  relation.  the  wave  scattering 
characteristics  of  the  location,  or  structural  junction,  can  be  actively 
altered  in  order  to  vary  the  path  of  power  transmission  or  reduce  power 
emanating  from  the  junction. 

Such  a  technique  has  several  advantages.  Wave  control  approaches 
the  problem  as  feedforward  disturbance  rejection  of  incoming  waves.  In 
some  applications,  this  can  eliminate  resonant  behavior  by  creating 
matched  terminations.  Local  models  are  insensitive  to  all  but  local 
modelling  errors.  Since  a  wave  model  is  a  local  description,  a 
guarantee  of  control  stability  is  not  based  upon  knowledge  of  global 
structural  behavior.  Instead.  stability  is  judged  based  on  the 
frequency  dependent  power  generation/dissipation  properties  of  the 
active  junction.  Global  control  performance  is  dependent  upon  the 
significance  of  the  disturbance  path  containing  the  active  junction. 
The  control  energy  expended  is  of  the  order  of  the  disturbance  energy 
and  only  a  few  actuators  and  sensors  are  required  to  carry  out  most 
tasks . 

Two  wave  control  formulations  have  been  developed.  The  first 
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involves  specifying  the  desired  scattering  behavior  and  deriving  the 
appropriate  control  The  second  minimizes  a  cost  based  on  the  power 
emanating  from  a  junction  and  the  level  of  control  effort  used.  These 
two  formulations  can  be  used  separately  or  in  combination  with  the 
specified  behavior  being  used  as  a  constraint  in  the  cost  minimization 
problem.  An  example  would  call  for  minimizing  power  emanating  from  an 
intersection  of  several  structural  members  while  prohibiting  waves  from 
departing  the  junction  along  one  of  the  members.  This  application  is  a 
combination  of  vibration  suppression  and  dynamic  isolation.  Both 
formulations  are  frequency  domain  techniques. 

When  employing  frequency  domain  techniques,  issues  of  causality 
must  be  dealt  with.  The  first  step  in  the  solution  of  the  cost 
minimization  problem  is  to  construct  an  integral  over  the  entire 
frequency  range  whose  integrand  is  composed  of  two  quadratic  terms.  The 
first  term  is  junction  power  flow  expressed  in  wave  mode  amplitudes  and 
the  second  term  is  control  effort.  Minimizing  the  integral  yields  an 
expression  relating  control  effort  to  wave  mode  amplitudes.  Junction 
deflections  can  be  substituted  for  wave  mode  amplitudes  to  provide  a 
causal  set  of  measurements.  The  resulting  relation  is  then  placed  in 
Wiener-Hopf  form  and  solved  accordingly  to  yield  the  optimal,  causal 
solution.  Alternatively.  some  researchers  have  proposed  methods  for 
estimating  wave  modes. 

For  some  applications,  such  as  the  matched  termination  of  a  rod. 
wave  control  resembles  direct  velocity  feedback.  In  general  the 
solutions  will  be  quite  different.  For  example,  solutions  for  beams  in 
transverse  bending  yield  feedback  of  half  derivatives  and  half  integrals 
of  junction  deflections  which  are  easily  implemen table  compensators  yet 
not  obvious  solutions  when  modal  models  are  employed.  In  theory,  these 
wave  controllers  yield  performance  far  exceeding  direct  velocity 
feedback  even  to  the  point  of  eliminating  resonant  behavior. 

Many  possibilities  stem  from  the  use  of  these  local  wave  models. 
The  compensators  may  be  adapted  based  upon  comparison  of  actual  outgoing 
wave  modes  with  their  predicted  levels.  Measurements  can  be  performed 
"upstream"  to  compensate  for  actuator  and  sensor  lags.  Vibrational 
energy  can  be  shunted  around  performance  critical  locations  or  trapped 
in  noncritical  regions  to  be  dissipated.  Many  of  these  possibilities 
are  not  readily  observed  using  modal  models  yet  are  obvious  using 
equally  valid  wave  models.  This  promises  to  be  cm  interesting  area  for 
research  which  addresses  many  of  the  performance  and  implementation 
requirements  necessary  for  large  space  structure  applications. 


REFERENCE 

[1]  Miller,  David  W. ,  von  Flotow,  Andreas  H. .  Hall,  Steven  R.,  "Active 
Modification  of  Wave  Reflection  and  Transmission  in  Flexible 
Structures,"  Proceedings  of  the  1987  American  Control  Conference,  Vol. 
2.  pp.  1318-1324. 


50 


Adaptive  Control  of  Large  Space  Structures 


Robert  L.  Kosut 
Integrated  Systems.  Inc. 
2500  Mission  College  Blvd. 
Santa  Clara,  CA  05051 


INTRODUCTION 

In  this  talk  we  briefly  describe  some  of  the  research  issues  involved  in  the  design  and  , analysis  of 
adaptive  control  systems  for  large  space  structures  (LSSi.  The  need  for  adaptive  control  arises  from 
many  envisioned  future  LSS  missions  which  impose  stringent  performance  demands  on  tracking 
accuracy  and  structural  vibration  attenuation.  Both  active  feedback  control  and  passive  damping 
will  be  i  practical  necessity,  and  moreover,  their  design  will  require  a  model  of  the  LSS  system  whose 
accuracy  is  compatible  with  the  performance  demands.  Structural  variations  from  many  sources, 
such  as  deployment,  material  fatigue,  and  even  random  variations  in  materials  and  manufacturing 
tolerances,  will  significantly  degrade  closed-loop  performance,  e.g  .  ACOSS(lOSl).  Thus,  the  on- 
orbit  dynamics  <. ;  LSS  will  not  be  sufficiently  like  those  obtained  from  either  ground-testing  or 
even  from  sophisticated  computer  generated  modeling  techniques,  such  as  finite  element  modeling. 
Current  structural  modeling  techniques  are  just  not  sufficiently  accurate  or  able  to  account  for  ail 
the  possible  sources  of  parameter  variation  Therefore,  under  these  conditions,  it  will  be  necessarv 
to  identify  the  LSS  dynamics  directly  from  on-orhit  measurements,  and  simultaneously,  tune  or  re¬ 
design  the  control.  Hence,  the  control  design  cycle  will  he  an  adaptive  process,  typically  starting 
with  a  nominal  has -performance  design  h.L-  d  on  a  coarse  model,  and  then  re-designed  from  on-orhit 
data.  The  adaptation  may  take  place  either  otf-hne  or  simultaneously  while  the  system  is  operating, 
but  in  either  case,  the  process  will  require  on-orbit  data. 


STRUCTURE  OF  ADAPTIVE  CONTROL  SYSTEM 

Adaptive  control,  as  depicted  in  figure  1.  consists  of  a  controller  and  two  additional  processes, 
namely:  (1)  a  model  estimator,  and  (21  a  control  design  rule  The  model  estimator  operates  on  the 
input-output  data  obtained  from  measurements  of  the  system  to  he  controlled,  producing  a  model 
estimate  from  a  pre-detennined  model  set  I  he  model  estimate  is  transformed  by  the  control  design 
rule  into  a  controller  to  he  used  in  feedback  with  the  actual  system  The  fundamental  question  is: 
when  will  it  work?  That  is.  under  what  conditions  will  the  model  estimate  converge  to  a  good  model 
in  the  allowed  model  set  By  a  good  model  is  meant  one  that  produces  a  controller,  via  the  control 
design  rule,  which  when  applied  to  the  actual  system  yields  acceptable  closed-loop  performance. 
Even  if  such  a  good  model  exists,  the  estimated  model  may  not  converge,  even  if  the  estimated 
model  is  initialized  close  to  a  good  model  Moreover,  if  convergence  is  too  slow,  then  unacceptable 
behavior  can  occur  during  the  learning  process.  1  bus.  convergence  alone  is  not  a  sufficient  condition 
for  establishing  good  performance  \Ve  also  need  information  about  t  he  convergence  rate,  particularly 
how  it  can  he  affected  by  user  choices  of  inputs,  data  filters,  etc. 

Ill  summary,  to  evaluate  an  adaptive  system,  it  is  necessary  to  determine  the  relations  among 
the  following  choices: 

(1)  a  representative  "true''  system,  (2)  a  parametric  model  set.  (3)  a  model  estimator  or  adaptive 
algorithm.  (1)  a  set  of  closed-loop  performance  criteria,  and  (5)  a  control  design  approach,  rule,  or 
algorithm. 
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Figure  1:  Adaptive  Control 


Mode!  Parametrization 

We  will  restrict  our  discussion  hero  to  parametric  transfer  function  models.  There-  are  two  basic 
choices  for  parameters  (1)  physical  parameters,  and  (2)  canonical  parameters.  Physical  parameters 
include  masses,  spring  constants,  electrical  circuit  quantities,  etc.  Canonical  parameters  are  best 
exemplified  by  transfer  function  coefficients,  which  usually  depend  in  a  complicated  manner  on  tie- 
physical  parameters  Moreover,  if  the  transfer  function  represents  a  sampled- data  system.  then  the 
relation  between  the  physical  parameters  and  the  transfer  function  coefficients  becomes  even  more 
complicated,  particularly  when  the  sampling  rate  is  low. 

In  parameter  estimation,  the  above  parametrizations  offer  different  adavantages  and  disadvan¬ 
tages.  The  estimation  of  transfer  function  coefficients  can  be  very  direct,  because  it  is  possible  to 
express  the  prediction  error  as  a  linear  function  of  these  parameters,  thus  leading  to  simple  estima¬ 
tion  algorithms,  such  as  recursive  least  squares.  Physical  parameters  are  directly  meaningful,  but 
because  they  enter  in  a  complicated  manner  into  the  transfer  function  coefficients,  the  prediction 
error  is  also  a  complicated  function  of  the  physical  parameters,  making  direct  estimation  schemes 
difficult  to  design.  Often  a  system  model  can  he  characterized  by  a  few  physical  parameters,  whereas 
the  transfer  function  may  consist  of  very  many  coefficients. 

One  can  now  ask  the  question:  which  of  the  above  mode]  paramelrizations  is  '‘better"?  To 
answer  this  question,  we  first  remark  that  ncif/ier  is  correct,  because  the  true  system  is  different 
than  either  model  Even  the  so  called  physical  transfer  function  model  is  undoubtedly  derived  under 
certain  assumptions,  such  as  small  deflections,  linear  elasticity,  uniform  mass  density,  etc.  Hence, 
t'-e  choice  should  be  made  by  considering  the  intended  use  of  the  model 

Model  Evaluation 

In  control  design,  which  is  the  case  we  consider  ln-re,  what  is  need-  d  is  an  estimate  of  the  transfer 
function,  whose  accuracy  is  determined  by  the  closed-loop  performance,  e  g.,  Doyle  and  Stein(lOSl), 
Safonov  et  al.fT'.WI).  Kosut  et  al,(I9S3),  Vidyasagar(1985).  Dhaya  and  IVseorf  19851.  Thus,  for 
cor, to!  design,  the  parametrization  can  be  viewed  merely  as  an  intermediate  step  However, 
different  parametric. items  invoke  parameter  estimation  algorithms  which  may  have  significantly 
dilT'-r-'nt  convergence  rates  and  numerical  stability  properites,  eg.  I.jung  and  Sodestromf  1983). 


Hence,  although  closed-loop  performance  is  the  ultimate  aim.  other  criteria  must  Lie  considered, 
such  a.s  those  pertaining  to  computation  of  the  estimate 

To  evaluate  the  model  paramet ri/ation  requires  the  f.-Howitig  fruits  from  I.jutld  1‘JST)  r<  warding 
parametric  transfer  function  estimation: 

•  Least  squares  parameter  estimation  is  equivalent,  asymptotically  in  the  length  of  the  data 
record,  to  minimization  of  a  weighted  frequency  domain  criterion  which  penalizes  the  quadratic 
error  between  the  true  system  transfer  function  anil  the  parametric  model  transfer  function 
The  weighting  function  depends  on  various  user  choices,  such  as  the  input  spectrum  data 
filters,  and  the  noise  or  disturbance  model. 

•  The  mean-square-error  (MSK1  between  tin  true  system  transfer  function  and  the  ideal*  es¬ 
timated  transfer  function  can  lie  explicitly  evaluated  and  depends  on  the  user  choices  listed 

above. 


One  of  the  most  significant  .aspects  of  the  above  results  is  that  the  MSB  litween  a  particular  choice 
of  model  paramet nzat ion  and  the  true  system  can  be  evaluated  witfunit  specifying  how  the  estimate 
is  computed  or  whether  i!  comt  ryes  I  iie  convergence  question  involves  the  analysts  techniques 
reported  in  Anderson  et  al  (  IfWi’O  IL-nce.  conditions  for  convergence  add  further  constraints  on 
the  user  choices,  it  is  important  to  emphasize  that  since  the  above  results  depend  on  the  choice 
of  a  "true''  representaive  system,  they  can  he  used  for  analysis  and  design  of  the  estimator  and 
estimation  experiment 

Many  research  efforts  have  been  undertaken  to  evaluate  methods  of  system  identification  for 
distributed  parameter  and  I.SS  type  systems,  e  g.,  to  name  a  few.  Goodson  and  Pol is(  197-1),  Banks. 
Crowley,  and  Kumscld  i'.'Mi).  Schaeciiteri  ltlsii).  and  Denman  et  al.(19Sl)l.  The  difference  betwe.-n 
these  reported  results  and  what  we  are  emphasizing  here,  is  the  connection  between  the  criteria  for 
system  identification  and  the  intended  use  of  the  estimated  mode]  for  control  design. 

As  one  example  of  the  kind  of  analysis  that  might  be  accomplished,  consider  the  use  of  simple 
continuum  models,  such  as  beams,  shells,  plates,  etc.,  with  parametrized  spatially  varying  param¬ 
eters.  Lor  example,  a  true  system  can  be  represented  as  a  non-uniform  elastic  rod.  whereas  the 
model  set  is  a  uniform  rod  with  the  parameters  to  be  estimated  being  the  constant  tortional  mass 
and  stiffness.  Iti  tins  case  it  is  possible  to  determine  what  amount  of  distributed  tortional  mass  and 
stiffness  can  be  sufficiently  well  approximmated  by  a  uniform  mass  and  stiffness  for  the  purposes  of 
control  design. 

The  above  simple  example  illustrates  the  difficulties  involved  in  choosing  a  model  parametrtza- 
tion.  Obviously  the  extent mn  to  the  BPS  environment  is  not.  so  straightforward  because  of  the 
additional  complexity. 

Experiment  Design 

The  MSB  can  be  used  to  aid  in  designing  the  estimation  experiment  when  the  intended  use  of  the 
estimated  model  is  control  design.  The  basic  idea  is  to  make  the  criterion  for  estimation  similar  to 
the  criterion  for  control  design  I  his  is  equivalent  to  minimizing  a  weighted  norm  of  the  MSP,  with 
respect  to  free  parameters  or  choices  jn  the  experiment,  such  as,  the  input  spectrum,  model  order, 
data  filters,  and  so  on,  as  desetbed  m  I,jung(  1 9S7) 

At  the  present  time  the  use  of  the  MSP  us  a  measure  of  experiment  design  for  t. fie  case  when 
tlie  intended  use  of  the  transfer  function  estimate  is  control  synthesis  is  still  in  the  beginning  stages 

1  I  hr’  ideal  estimate  ininimr/.es  the  least  Mjures  rnternm  for  a  specified  length  of  the  data  record.  This  is  not  the 
estimate  which  is  recursively  computed  fi  om  on  line  data.  These  two  estimates  are,  at  best,  asymptotically  equivalent, 
with  the  rei  ursive  estimate  having  a  larger  variance,  the  si/.e  of  which  is  dependent  on  the  specific  algorithm  for  on-line 
estimation. 
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One  of  the  major  problems  is  that  measures  of  closed  loop  performance  are  complicated  functions 
of  the  MSI-,  and  the  criterion  for  estimation 


Computing  Model  Uncertainty  from  Data 

The  requisite  information  for  robust  feedback  stabilization  is  typically  a  nominal  model  of  the  plant 
together  with  an  uncertainty  profile.  From  the  above  discussion,  computation  of  the  MSF  for  model 
error  allows  for  art  analyst*  of  the  adaptive  system.  If  the  MSE  can  be  estimated  from  on-line  data, 
then  a  robust  controller  can  be  designed  Our  current  approach  to  estimating  model  uncertainty 
involves  not  only  parametric  estimation,  hut  also  non-parametric  (spectral)  methods  of  transfer 
function  estimation,  e.g.,  Jenkins  and  \Vatts(  1968).  The  non-parametric  methods  provide  estimat*  s 
of  model  accuracy  as  a  function  of  frequency,  particularly  over  those  frequency  ranges  wlu-re  either 
the  model  structure  is  poorly  known,  or  else  the  model  accuracy  in  that  range  is  unimportant  for 
control  design.  Hie  feasibility  of  using  spectral  techniques  for  the  estimation  of  model  error  can  he 
found  in  Kosut(  1087). 


Stability  and  Convergence  Analysis 

Conditions  for  stability  and  convergence  and  causes  for  divergence  and  instability  in  adaptive  control 
systems  are  contained  in  Anderson  et  al.(ldSC)  The  stability  analysis  focuses  on  slow  adaptation 
and  the  use  of  some  of  the  classical  methods  for  analyzing  differential  equations,  e  g.,  linearization, 
the  method  of  averaging,  and  Lyapunov's  second  method,  see,  e.g.,  Hnle(1969).  The  material  in  the 
text  represents  modifications  and  refinements  of  earlier  work  by  the  authors. 

To  remove  the  restrictiveness  of  slow  adaptation  requires  an  understanding  of  the  transient 
behavior  of  adaptive  systems.  Preliminary  investigations  are  reported  in  Kosut  et  aid  1987), 
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INTRODUCTION  AND  BACKGROUND 

The  problem  addressed  here  is  the  determination  it  bounds  in  the 
degradation  of  system  performance  due  to  uncertainties  and/or  ur.’a  r<  a  ter. 
and  it. perfectly  modelled  subsystem  interact  ions.  Such  bounding 
techniques  represent  a  fundamental  systems  «ir.  >i  >  sis  tool  that  is 
indispensable  for  further  elucidation  ot  decentralised  controller 

ci  L  C  li  Z.  t.  <L‘C  C  U  r  fc*£>  ciliu  1  C'DU  J  t 

Extensive  word  has  been  carried  out  within  the  controls  community 
ir.  the  urea  of  f :  t-q  none  y  -  d  cm  a  in  analysis  ol  robust  .  t.ibility  1  iv  it.,,  vise 
to  th.e  H-infir.ity  theory  of  robustness  characterization  and  robust 
design  [1-5].  However,  on  several  occasions  we  l.uvt  remarked  that 
although  the  K-infinity  world-view  is  a  beautiful  and  compelling  theory 
within  its  proper  province,  its  fundamental  assumptions  lender  it 
inapplicable  to  structural  vibration  control  which  involves  parametric 
and  often  nondestabilizing  open-loop  uncertainties.  A  principal 
difficulty  is  the  conservatism  of  H-infir.ity  robustness 

characterizations.  A  stability  robustness  analysis  technique  is  called 
conservative  if  the  predicted  set  of  nondestab  il  izir.g  perturbations  is  a 
proper  subset  of  the  ac  tual  set  of  r.ondestabil  izing  perturbations.  Note 
that  conservatism  jointly  depends  upon  both  the  definition  cf  admissible 
perturbation  classes  and  the  robustness  analysis  technique. 

The  well-known  conservatism  of  H-ir.finity  theory  does  not  arise 
because  it  operates  in  the  frequency  domain,  per  se,  or  because  the 
infinity  norm  is  employed,  but  rather  because  of  the  crudeness  cf  K- 
infinity  bounds.  What  is  the  fundamental  source  or  this  crudeness? 
Possibly  this  arises  because  the  fundamental  intent  of  K-infinity 

development  was  the  extension  of  classical  control  design  concerts _ -v 

the  mu ltivariable  case  whether  or  not  classical  concepts  c.re  truly 
suited  to  the  problem  at  hand. 

For  example,  in  keeping  with  classical  ideas,  there  has  beer 
widespread  insistence  upon  couching  all  questions  or  performance  and 
uncertainty  in  terms  of  simplistic  (albeit  traditional)  unity  gain 
feedback  diagrams.  Thus,  singular  value  developments  have  lumped 
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uncertainty  in  a  single  block  thereby  obscuring  the  often  complex 
structure  of  modelling  error.  Moreover,  this  feedback  paradigm  is 
maintained  even  for  structured  uncertainty  approaches  [6]  . 

LARGE  SCALE  SYSTEM  FORMULATION  AND  Lp  BOUNDS 

Here,  we  contend  that  to  achieve  a  less  confining  point  of  view 
the  tirst  stop  is  to  represent  uncertain  systems  by  means  of  a  large- 
scale  system  input-output  formulation  as  depicted  in  Figure  1. 
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Figure  1.  Large-Scale  System  Input -Output  Formulation 
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Referring  to  Figure  1,  the  overall  system  is  represented  by 
interconnected  subsystems  undergoing  interactions.  The  subsystems, 
characterized  by  the  operators  G,  (k=l . r) ,  represent  the  known 

dynamics  of  the  system  while  the  subsystem  interactions,  given  by  the 

operators  H.  correspond  to  uncertainties.  Note  that  tne  partitioned 
kj 

off-block-diagonal  operator  H  is  stipulated  to  belong  to  some  compact 
arcwise  connected  set  H.  The  set  H  specifies  both  the  character  and 
extent  of  dynamical  uncertainties. 

The  motivation  for  the  above  input-output  formulation  within  the 
context  of  large-scale  systems  is  obvious.  But  in  ••edition,  thanks  to 
the  Dynamic  Inclusion  Principle  and  related  ideas  elaborated  by  Siljak 
and  his  co-workers  [7,8]  the  representation  of  Figure  1  is  also  suitable 
for  parametric  perturbations  in  monolithic  systems,  i.e.,  systems 
without  explicit  interconnections. 

The  p rcblem  now  addressed  is  how  to  bound  the  degradation  of  the 
system  output  vector  y  or  the  prediction  accuracy  y-y  due  to  the 

uncertainties . 

To  give  this  problem  mathematical  form,  we  must  use  the  block- 
matrix  results  of  Ostrowski  [9]  and  define  the  block-Lp  norm  matrix:  of  a 
partitioned  operator,  as  the  nonnegative  matrix  whose  elements  are  the 
Lp  norms  of  the  corresponding  subblocks.  With  this  definition,  the 
principal  problem  is  to  bound  the  block -rortr.  vector  of  the  system  output 
y  ever  all  variations  of  the  uncertain  perturbations.  Note  also  that 
the  existence  of  such  a  finite  bound  implies  input-output  stability 
(see  [10] ) . 

Now,  the  above  articulation  of  uncertainties  into  numerous 
interactions  permits  more  finely  articulated  methods  of  computing  bounds 
beyond  singular  value  analysis,  namely,  methods  associated  with  the 
majorant  analysis  of  Dahlquist  [3.11]. 

The  uncertain  subsystem  interaction  format  (introduced  in  Figure  1] 
in  conjunction  with  majorant  analysis  gives  an  almost  unlimited 
potential  for  formulating  sharper  bounds.  Using  a  process  of  operator 
iteration,  one  can  obtain  a  hierarc hy  of  output  bounds,  where  each 
successive  member  of  the  hierarchy  requires  more  and  more  information 
but  is  less  and  less  conservative  (with  respect  to  the  set  H) .  Also, 
because  we  work  in  an  operator  setting,  distinctions  between  the  time 
and  frequency  domains  are  blurred.  It  is  parochial  to  assert  that  only 
frequency-domain  or  time-domain  methods  must  be  used.  What's  needed  is 
easy  and  fluent  translation  between  the  frequency  and  time-domain. 
Furthermore,  the  computational  advantage  of  this  kind  of  hierarchy  is 
that  each  bound  requires  only  the  inversion  of  an  M-matrix.  This  is 
quite  straightforward  and  nicely  tractable,  even  for  many  subsystems, 
since  it  involves  computing  a  monotonically  increasing  sequence  where 
each  iteration  involves  an  addition  and  a  multiplication  of  low-order 
matrices. 
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BOUNDS  FOR  RESPONSE  TO  STOCHASTIC  INPUTS 


The  above  discussion  has  sketched  the  general  development  of 
majorant  robustness  analysis  within  an  operator  setting  which  employs  Lp 
norms  to  describe  the  "size"  of  subsystem  outputs.  For  systems  with 
stochastic  inputs  and  time  independent  parameter  uncertainties,  the  main 
lines  of  development  are  analogous.  However,  in  this  case  one  needs  to 
work  with  the  Lyapunov  equation  for  the  steady-state  second-moment 
matrix  of  response  and  then  derive  majorant  bounds  for  the  block-norm 
matrix  or  the  second  moment.  The  general  setup  for  undertaxing  majorant 
analysis  for  parametrically  uncertain  stochastic  systems  is  shown  in 
Figure  2.  here,  the  block-diagonal  matrix.  A  represents  the  known 
subsystem  or  nominal  system  dynamics  while  the  of f -block-diagonal  matrix 
G  represents  uncertain  subsystem  interactions  or  parametric 
uncertainties . 
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Figure  2.  Subsystem  Interaction  Model 


Note  that  the  disturbance  intensity  V  and  the  second-moment  matrix 
Q  are  partitioned  conformably  with  A  and  G.  We  bound  performance 
degradation  due  to  uncertain  interactions  G  ranging  over  the  admissible 
set  G  by  bounding  the  block-norm  matrix  of  Q.  To  do  this,  however, 
requires  additional  algebraic  tools,  such  as  the  Block-Kronecker  Algebra 
which  is  a  generalization  of  the  matrix  calculus  reviewed  by  Brewer 
[12].  Because  of  the  algebraic  complexity  of  deriving  majorants  iur  the 
second-moment  matrix,  the  Krcnecker  Algebra  is  far  more  than  a  mere 
notational  convenience. 

In  particular,  using  the  Block-Kronecker  Algebra  one  first  reduces 
the  problem  to  a  more  tractable  form  and  then  applies  majorant  analysis 
to  obtain  a  hierarchy  of  majorant  bounds.  As  in  the  Lp  bound  analysis, 
each  successive  member  of  the  hierarchy  offers  less  and  less 
conservative  bounds. 
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We  now  consider  in  more  detail  the  first  two  members  of  the 
majorant  hierarchy  in  order  to  illustrate  the  specific  forms  of  the 
modified  Lyapunov  equations  that  are  obtained. 

For  example.  Figure  3  shows  the  first  member  of  the  second-moment 
majorant  hierarchy  [13]  .  This  gives  the  majorant  Q  as  the  solution  of  a 
simple  nonnegative  matrix  equation,  where  *  denotes  the  Hadatr.ard 
(element  by  element)  product  ar.d  the  row  and  column  dimension  of  the 
equation  is  the  number  of  subsystems.  For  the  norm-bounded  uncertainty 
set  shown  in  Figure  3,  the  existence  of  a  nonnegative  solution  implies  a 
bound  for  the  block-norm  matrix  of  the  second  moment  and  robust 
stability,  : . e . ,  (A+G)  is  stable  for  all  perturbations  G  in  the  norm- 
bounded  set. 


£ 


■  Robust  Stability 

■  Robust  Performance 

Figure  3.  Majorant  Lyapunov  Equation 


One  particular  advantage  of  the  first  member  cf  the  hierarchy  is 
that  it  correctly  shows  the  effect  of  wide  frequency  separation  of 
subsystems  on  performance  degradation  and  robust  stability.  In 
particular,  the  majorant  equation  will  correctly  predict  that  as 
frequency  separation  becomes  sufficiently  large,  subsystems  become 
effectively  decoupled.  Such  predictions  cannot  be  made  by  the  small 
gain  theorem  for  large-scale  systems  or  by  vector  Lyapunov  theory  [7,8]. 
Thus,  even  the  first  member  cf  the  majorant  hierarchy  offers  greatly 
reduced  conservatism  compared  to  previous  results. 

Furthermore,  the  second  member  of  the  second-moment  majorant 
hierarchy,  shown  in  Figure  A  gives  even  tighter  bounds  and  can  even 
predict  the  stabil i zing  effect  of  certain  kinds  of  perturbations.  The 
form  of  the  majorant  equation  (top  of  Figure  A)  is  similar  to  the  first 
member  of  the  hierarchy  except  that  the  operation  H [0]  appears.  This 
operator  is  precisely  what  would  arise  in  the  equation  for  the  second- 
moment  matrix  for  a  system  with  Stratonovich  noise  parameters!  Fo  far, 
we  have  discussed  a  design  analysis  tool  for  predicting  performance 
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degradation  due  to  uncertainty.  This  crucial  observation  brings  us  to 
consideration  of  the  link  between  majorant  robustness  analysis  and  MEOP 
design  synthesis  theory. 


Second  member  of  the  hierarchy: 
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■  Nominal  '  performance,  tr  [QR],  given  by  Stratonovich  model 

Figure  4.  Second  Member  of  the  Majorant  Hierarchy 

THE  LINK  BETWEEN  ANALYSIS  AND  SYNTHESIS 

Figure  5  illustrates  this  link  and  the  accompanying  sequence  of 
logical  developments.  Overall,  one  ray  regard  the  MEOP  design  synthesir 
theory  as  arising  from  a  particular  robustness  analysis  tool.  Although 
any  member  of  the  second-moment  majorant  hierarchy  might  be  chosen  as 
the  basis  of  a  design  synthesis  theory,  we  choose  the  second  member  of 
the  hierarchy  (see  lower  right  block  in  Figure  5)  to  serve  as  the  point 
of  departure  because  it  is  the  simplest  bound  that  also  handles 
nondestabilizing  uncertainties.  Referring  to  the  lower  left  block  of 
Figure  5,  it  is  seen  that  the  second-moment  equation  of  a  multiplicative 
Stratonovich  noise  model  essentially  gives  an  approximation  to  the 
majorant  equation  and  a  smooth  optimization  problem.  The  Stratonovich 
second  moment  equation  then  leads  to  an  auxiliary  optim ization  problem 
(upper  left  block  in  Figure  5),  namely,  choose  dynamic  compensator  gains 
to  minimize  the  quadratic  performance  of  a  system  having  mul t ipl icat ive 
stochastic  parameters.  Because  of  the  Stratonovich  modifications  to  the 
standard  form  of  the  Lyapunov  equation  that  appear  in  the  equation  for 

Q,  the  robust  stability  condition  implied  by  the  majorant  equation  is 
still  enforced  since  the  optimization  problem  imposes  a  robust 
performance  constraint. 


Auxilliary  Optimization  Problem 
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Figure  5.  Majorant  Hierarchy  and  Stratonovich  Models 
—  Link  Eetween  Analysis  and  Synthesis 

This  optimization  of  an  apparently  stochastic  system  actually 
approximates  the  majorant  bound  which  was  derived  purely 
deterministically  and  leads  to  the  rather  elegant  fiEOP  optimality 
conditions  given  in  the  upper  right  block  in  Figure  5. 

Of  course,  the  use  of  Stratonovich  stochastic  models  was  earlier 
indicated  by  maximum  entropy  principles  and  stochastic  approximation 
theory,  and  this  line  of  development  still  stands.  But  the  import  of 
the  more  recent  majorant  analysis  developments  is  that  there  is  a  direct 
1 mk  between  maximum  entropy  stochastic  modelling  and  deterministic 
performance  bounds.  This  link  tends  to  blur  the  distinctions  between 
stochastic  and  deterministic  points  of  view.  This  is  just  as  well:  The 
task  confronting  the  controls  and  systems  theory  community  is  not  to 
resolve  the  stochastic  versus  deterministic  debate  one  way  or  the  other, 
but  rather  to  rise  above  it.  As  the  work  described  here  suggests,  there 
is  a  plane  upon  which  these  points  of  view  are  numerically 
indistinguishable . 
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INTRODUCTION 

Despite  significant  advances  in  the  cost  ana  pc-rtcrmunce  it  digital 
computers  ever  the  last  decade,  there  remains  a  need  in  several 
technological  areas  for  low-order,  high-per f ornate e  controllers,  Ir. 
particular,  this  paper  is  motivated  by  the  problem  of  vibration 
suppression  in  large  flexible  space  structures.  Such  systems  are 
i  r.  f  in  ite-dimenuicnal  (distributed  parameter)  in  nature  and  her.ee  <»ny 
finite-dimensional  controller  is  necessarily  of  reduced  order.  The  need 
for  low-order  controllers  is  further  driven  by  severe  constraints  on 
cost,  weight  and  power  in  space  systems,  not  to  mention  the  restriction 
to  space-qualified  computational  hardware. 

A  wide  variety  of  approaches  have  been  proposed  to  obtaining 
reduced-order  controllers.  A  comparison  of  several  approaches  to 
controller  reduction  is  given  in  [1],  These  metho  ir  operate  by  first 
designing  a  high-order  LQG  controller  and  then  obtaining  a  suitable  low- 
order  controller  by  means  of  controller  reduction. 

A  more  direct  approach  to  designing  reduced-order  controllers 
involves  optimizing  the  quadratic  performance  functional  over  the  class 
of  controllers  of  fixed  order.  The  controller  order  may  be  determined 
by  implementation  constraints  or  can  be  varied  for 
per formanc e/ throughput  tradeoff  studies. 

A  reformulation  of  the  parameter  optimization  approach  was  given 
recently  in  [2].  The  authors  showed  that  the  first  order  necessary 
conditions  can  be  transformed  to  yield  expl  ic  i  t  gain  t-xj  rest  ions  lor 
extremal  fixed-order  controllers. 

Fegardless  of  how  appealing  the  optimal  prcjecti.r.  :  :t.ulut  iij. 
[2-12]  may  appear  to  be,  its  contribution  is  vacua:;  in]  ess  certain 
serious  questions  car.  be  resolved.  Those  include: 

1.  Under  what  conditions  on  the  problem  data  car.  the  optimal 
projection  equations  be  guaranteed  a  priori  to  possess  a 
solution? 
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2.  Given  problem  data,  exact Iv  how  many  solutions  do  the  equations 
possess? 

3.  Or  the  possible  solutions,  wnat  are  their  stability  properties, 
what  is  their  performance,  and  which  is  the  global  optimum? 

4.  How  can  numerical  algorithms  be  constructed  which  can  be 
guaranteed  to  converge  to  any  desired  solution  especially  the 
global  minimum? 

It  seems  clear  that  any  attempt  to  address  the  above  issues  aiust 
utilise  mathematical  methods  which  are  global  in  nature.  To  this  end  we 
have  applied  degree  theory  and  associated  homotopic  continuation  methods 
([13-24])  to  analyse  the  solutions  to  the  optimal  projection  equations 
and  to  construct  convergent,  implo-ner.table  algorithms  for  their 
coniDu ta t i on .  The  purpose  of  this  presentation  is  to  report  significant 
recent  results  in  this  regard. 

HOMOTOPIC  CONTINUATION  AND  DEGREE  THEORY 
FOR  THE  OPTIMAL  PROJECTION  EQUATIONS 

A  homotopic  continuation  method  for  solving  a  problem  is  to  first 
solve  an  tasy  "similar"  problem,  uni  then  to  continuously  deform  the 
easy  problem  into  the  original  problem  and  to  follow  the  path  of 
solutions  as  the  easy  problem  is  deformed  into  the  original  problem. 

This  is  shown  conceptually  in  Figure  1. 
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Figure  1 


Topological  degree  theory  is  the  basic  tool  tor  analysing  the 
existence  and  topology  of  continuation  solutions. 


Here  we  consider  how  these  concepts  have  been  recently  applied  by 
Richter  [25]  to  solution  of  the  Optical  Projection  Equations  (OPE).  The 

basic  computational  problem  is  to  find  P,  Q,  P,  Q  which  satisfy: 
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For  t=0  the  solution  is  easy  to  find.  The  object  is  to  follow  the  path 
or  paths  of  solutions  F(t),  Q(t),  P(t),  Q(t)  from  t  =  0  to  t=l. 

In  following,  these  initial  solutions  from  t=0  to  t  =  l  there  are 
several  situations  which  could  occur  (see  figure  2). 

Using  degree  theory,  Richter  [25],  has  show’n  the  following  results: 

Theorem  1 


With  respect  to  the  above  continuation  method  of  OPE  solution,  the 
situations  shown  in  dashed  lines  in  Figure  2  cannot  occur.  That  is,  the 
only  solutions  to  the  OPE  at  t=l  (or  for  C^t<_l)  are  those  which  are 
continuously  derived  from  the  solutions  at  t=0. 
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t  =  u  i  1  =  1 


Topological  degree  theory  The  dashed  paths  do  not  exist 

If  Nc  >  mm  >\  N)-n.  there  is  only  one  solution  (nj  global  minimum) 


Theorem  2 


F i gure  2 


Assume  that  the  plant  is  stabilizable  and  detectable,  >  C, 

R.  >  0  and  r:  <  n  when  n  is  the  dimension  of  the  unstable  subspace  of 
1  u  -  c  u 

A.  Then,  in  the  class  of  nonr.egative-def inite  solutions  0,  P,  Q,  P  with 


rank  Q  =  rank  P  =  rank  QF  -  n  , 

c 


the  optimal  projection  equations  possess  at  most 


(min(n,m,i)  -n 

U 

r.  -n 
c  u 


) 


n  <  min  'n.tr.A)  , 
c  - 


1 


otherwise. 


stabilizing  solutions.  Each  such  solution  is  reachable  via  a  homotopic 
path  with  starting  point  corresponding  to  diagonal  initial  cata. 
Furthermore,  if  the  plant  is  stabilizable  by  means  of  an  n^th-order 

dynamic  compensator,  then  there  exists  at  least  one  solution. 

Note  that  when  n^  is  larger  than  the  number  of  inputs  or  outputs, 

there  is  only  one  solution  to  the  CPE  and  this  solution  corresponds  to 
the  global  minimum  ot  the  quadratic  performance  index. 
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ALGORITHM  DESCRIPTION  AND  NUMERICAL  RESULTS 


In  a  homotcpy  path  following  algorithm  one  follows  the  path  of 
solutions  oi  F(x(t),t)-0  by  integrating  the  initial  value  problem 

r’F 

-j-  =  F  (x  ( t ) )  j- (x  ( t ) ,  t)  ;  x(0)  =  x  . 

ut  x  dt  0 

For  the  optimal  projection  equations  the  solution  P.  P,  Q,  Q  can 

be  easily  determined  once  the  projection  r  is  known  so  that 

?(t)  -  F  ( r  ■  t ) ) ,  Q(t)  =  Q  ( -^  ( t ) )  etc.  Titus  the  derivatives  or  P,  Q,  P,  Q 

can  be  written  in  teres  of  derivatives  of  G‘  and  /"where  G1  and  P  ..re 

the  ractors  of  r{r  -  Gp)  .  Thus  we  obtain 


which  gives  2n  n  equation  for  P'  and  G*  .  F',  Q',  ?'  and  Q'  .-.re  ti  er. 

T*  i 

calculated  from  r'  and  G*  and  finally  p t  +4t )  is  updated  by 

P\ t+4t)  =  Tit)  +  P'  x  4t 


and  likewise  for  G,  P,  Q,  P  and  Q. 


Mote  that  each  continuation  step  involves  the  solution  of  a  1  inear 

system  ot  dimension  In  n.  This  means  that  for  small  controller  din.en- 

c 

sicn  ii.  the  c cciputat ional  burden  of  the  OPE  homotopy  algorithm  can 

actually  be  less  than  the  LQG  reduction  methods  (which  typically  require 
two  nxn  F.iccati  and  two  nxn  Lyapunov  equation  solutions). 


Figure  3  summarizes  the  closed-loop  stability  results  reportc-d  in 
[1]  for  LQG  reduction  methods  along  with  results  obtained  using  the 
homotopy  method  for  solving  the  optimal  projection  equations.  Here  qn 

is  a  scale  factor  for  the  plant  disturbance  noise  affecting  controller 
authority.  Clearly,  LQG  reduction  methods  experience  increasing 
difficulty  as  authority  increases,  i.e.,  as  the  r ^  terms  become 

increasingly  tore  important  m  coupling  the  control  and  reduction 
operations.  Ir.  contrast,  the  OPE  consistently  yield  closed-loop  stable 
reduced-order  controllers  for  all  values  of  q„.  Figure  4  shows  the 

quadratic  performance  ,J,  versus  the  observer  bandwidth  parameter  q„  ror 


OPE  designs  of  various  ordt'rs.  Ever,  lot  very  high  control  authority 
(q^-lO  ).  when  LQG  reduction  methods  produce  no  stabilizing  designs,  the 
lid 

2  order  CPE  compensator  shows  a  performance  only  20%  higher  than  a 
full-order  LQG  design. 
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File  complexity  ol  present  and  envisaged  space  structures  dictates  t  lie  inevitable 
need  lor  emphasis  on  cont  rol  and  strucliire  interact  ion.  More<  >ver.  in  view  of  present 
methodologies,  stringent  control  requirements  such  as  precise  pointing  and  slewing, 
vibration  suppression,  .and  shape  control  indicate  1 1 1 a  1  much  work  remains  to  be 
done  in  this  important  area.  In  this  presentation  we  summarize  recent  results  in 
tin’  decentralized,  relegated  control  oi  huge  space  structures,  with  concentration  on 
topics  ol  decentralization,  relegation,  servomechanism  design,  and  multiple  mirror 
system  examples. 

Decentralization  and  Relegation 


(Jeneral  goals  of  this  research  have  been  focused  in  the  area  of  large  space  structure 
(  LSS  )  modeling  and  cont  rol.  In  particular,  t  he  work  Inis  been  orient  eel  to  evaiuat  ing 
the  use  oi  decentralized  control  strategies  as  a  fundamental  context  within  which  to 
develop  controller  design  methodologies  for  a  varietv  of  control  objectives  for  large 
flexible  struct  tires.  It  is  within  this  context  that  we  develop  decent  ralized  servocom- 
pensator  design  techniques  which  incorporate  various  implementation  constraints, 
such  as  subsystem  input  bandwidth,  sensitivity  and  robustness  requirements,  in¬ 
tegrity  constraints,  and  so  on.  The  concept  of  relegation  as  coordinated  with  de¬ 
centralized  schemes  is  also  being  developed.  Relegation  is  the  assignment  of  control 
tasks  and  information  channels  in  view  of  control  efiecti veness.  on-line  computa¬ 
tional  complexity,  controller  capabilities,  and  plnsical  and  structural  constraints. 
Based  on  the  system  model  formulation  and  proposed  control  objectives,  we  believe 
relegation  to  be  ext  reiiie!  y  import  ant  for  the  cont  rol  of  large  scale  s\st  ems  in  general 
and  of  large  space  structures  in  particular. 
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Advanced  space-based  detection,  surveillance,  and  defense  systems  necessitate  the 
const  ruction  of  LSS  with  hip;  hi  v  flexible  component  s.  Such  s\  st  ems  must  opera  t  e  ef¬ 
fectively  under  disturbances  (for  example,  constant  cli  st  urbances  due  to  solar  winds 
or  sinusoidal  disturbances  due  to  the  operation  of  equipment  on  board  )  and  must  be 
able  to  track  moving  targets  or  reference  trajectories.  Difficulties  arise  from  the  lack 
of  significant  inherent  damping  characteristics  and  from  the  fact  that  it  is  usually 
not  possible  to  obtain  precise  linear-time-invariant  models  lor  such  systems.  I  here- 
fore.  it  is  necessary  to  desipn  on-line  controllers  which  are  robust  to  uncertainties 
in  modeling,  to  achieve  desired  damping,  disturbance  rejection  and  tracking. 

In  the  case  of  I,SS.  the  solution  of  the  servomechanism  problem  may  require  decen¬ 
tralized  controllers  for  various  reasons  such  as: 

a)  1  he  sensor  actuator  locations  on  the  structure  may  have  a  decentralized  tia 
tin  •  which  makes  it  hard  or  impossible  to  implement  feedback  through  a 
centralized  location. 

b)  Information  transfer  may  be  costly  so  that  the  designer  mav  wish  to  impose 
decent  ralized  feedback. 

e)  On-line  centralized  feedback  calculations  may  be  too  time  consuming  so  that 
simple  local  feedback  calculations  may  be  preferred. 

A  LSS  can  in  general  be  described  by  a  generalized  wave  equation  which  leads  to  an 
infinite  dimensional  st  at  e-space  model.  Alt  ernat  i  vely.  finite-element  met  hods  can  be 
employed  to  obtain  a  finite  dimensional  model  direct |y.  However,  an  accurate  model 
may  necessitate  a  large  number  of  states,  which  makes  it  very  hard  or  impossible  to 
design  a  controller  based  on  such  a  model.  In  this  case,  it  is  necessary  to  use  model 
reduction  techniques  to  obtain  a  model  of  manageable  dimension. 

Decentralized  LQG  LTR 

I'lie  above  controller  design  problem,  which  is  usually  referred  to  as  the  robust 
servomechanism  problem  has  been  studied  by  various  researchers  jn  a  general  con¬ 
text.  The  linear-quadratic-Cau.ssian  with  loop-transfer-recovery  I  LQG  LTR)  de¬ 
sign  methodology  has  been  introduced  recently  and  has  attracted  much  interest. 
This  methodology  enables  frequency  domain  controller  design  bv  using  the  well 
known  LQG  techniques. 

Recently,  a  linear-quadratic  optimal  approach  has  been  developed  to  solve  this 
problem,  win  re  it  has  been  shown  that  the  resulting  font  roller  n  jovs  manv  desirable 
properties.  In  this  presentation  we  consider  this  development  with  the  LQG. 'LTR 
methodology  for  decentralized  control  systems.  Specificall v.  bounds  on  the  norm 
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of  the  multiplicative  error  matrix  are  calculated  ami  the  1.0(1  I.i  K  met  hodology  is 
used  to  design  local  controllers. 

Vibration  Compensation  in  Optical  Tracking  Systems 

The  focus  of  this  work  is  on  compensation  for  flexibility  elici  ts  due  to.  lor  example, 
slew  maneuvers,  for  flexible  structure  systems  composed  <>(  multiple  actuated  mir¬ 
rors.  Primary  applications  for  such  systems  include  line  <4  sight  pointing  systems 
on  large  flexible  structures,  space  telerobot ic  systems,  and  space  telescope  -systems. 

I  wo  stages  characterize  the  modeling  problem  tor  this  -tndv:  1  )  I  he  description  of 
the  rigid  slewing  motion  and  associated  mirror  and  ray  optics;  and  I  I  lie  descrip¬ 
tion  of  the  fieribh  dynamics.  The  rigid-motion  ray-trace  equations  are  developed 
by  using  the  compact  notation  used  to  describe  the  motion  of  robotic  manipulators, 
while  flexible  dynamics  are  obtained  via  standard  finite-element  techniques.  The 
resulting  hybrid  model  is  suitable  tor  analysis  in  a  four  stage  process:  1  )  Relegation 
of  control  tasks  (intimately  related  to  the  kinematicsl:  I  The  standard  slewing 
control  problem:  -I)  Flexibility  compensation  using  mirror  actuators;  and  1)  Active 
vibration  damping  with  additional  t  proof  mass|  actuation.  In  this  presentation  we 
describe  results  addressing  the  first  and  third  .'tapes  of  the  above  process. 
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INTRODUCTION 

There  are  many  obvious  reasons  why  simple  controllers  of  low  order  are  to  be  pre¬ 
ferred  to  high  order  controllers.  However,  established  controller  synthesis  techniques,  such 
as  LQG.  frequency  weighted  LQG.  and  H^  minimization  result  in  controllers  of  high  order, 
even  of  order  higher  than  the  plant.  While  this  may  be  tolerated  for  systems  of  small 
dimension,  it  is  obviously  unacceptable  where  the  plant  itself  is  of  high  order,  such  as  it  is 
with  even  the  simplest  space  structures. 

There  is.  therefore,  the  desire  to  design  low-order  controllers  for  high  order  plants. 
Procedures  to  solve  the  low-order  controller  design  problem  can  broadly  be  divided  into 
two  classes  [  1  ]:  Direct  methods,  in  which  the  parameters  defining  a  low  order  controller 
are  computed  by  some  optimization  or  other  procedure,  and  indirect  methods,  in  which  a 
high  order  controller  is  found  hrst.  and  then  a  procedure  is  used  to  simplify  it.  Lxamples 
of  direct  methods  are  the  parametric  LQ  designs  [5]  and  the  projective  controls  procedure 
developed  by  the  authors  [2.3.4].  In  the  case  of  indirect  methods  the  LQG,  the  frequency 
weighted  LQG.  and  the  H^-norm  minimization  approaches  provide  the  high  order  controller 
and.  each  in  its  own  right,  captures  many  relevant  performance  robustness  design  goals. 
However,  as  stressed  in  [l],  procedures  for  reducing  high  order  controllers  to  low  order  con¬ 
trollers  are  less  well  developed.  In  [l]  the  controller  reduction  problem  is  reduced  in  princi¬ 
ple  to  the  constrained  L^  model  matching  problem: 

K*r(s)  =  arg  min  li  Tj(s)— T,(s)  K(s)  T-,(s)  H  rj0  U  ) 

Kr(s)  of  order  r  or  less 

and  Kr(s)  is  the  desired  reduced  order  controller.  Unfortunately,  there  is  no  convenient 
algorithm  for  the  solution  of  the  constrained  model  matching  problem  in  the  or  L„ 
spaces. 


REDUCTIONS  AND  L^- BOUNDS 

Important  relations  have,  however,  been  established  between  certain  state  space  model 
reduction  techniques  and  optimal  H,n  and  L^,  approximations  that  are  useful  in  the  solution 
of  the  controller  reduction  problem.  They  will  be  exploited  here  in  a  novel  approach  which 
can  be  used  in  direct  methods  as  well  as  in  indirect  methods  of  solving  the  low-order  con¬ 
troller  design  problem,  and  in  many  other  control  problems.  For  the  controller  reduction 
problem  these  results  will  essentially  allow  one  to  minimize  a  bound  y.  and  to  ascertain 
that  II  T,(s)— T2(s)K(s)T3(s)  II  ^y. 

To  highlight  the  ideas  let  A.B.C  be  a  balanced  realization  [6]  of  a  linear  time-invariant 
strictly  proper  system  with  transfer  function  G(s).  in  the  sense  that 

G(s)  =  C(sl— A)1  B  (2) 


AE  +  IAT  +  BBl  =  0 


(3) 
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(4) 


AT£  +  LA  +  cTc  =  a 

where  I  is  the  (balanced)  observability  and  controllability  grammian.  with 

I  =  diagicr j.cr j . <xn|  .  cr  >  tr1+1  .  i  =  1 . n— 1 

It  is  known  that  the  k-dimensional  reduced  order  model 

G,(s)  =  C,(sl— A,  j)_1Bj 
where  Ajj.Bj.C,  come  from  the  partition 


•^n 

At]  A  22 


.  C  =  [C,  Ct] 


enjoys  the  following  properties: 

(i)  Svstems  (A  .B  .C  ).  i=l  .2  are  asvmptoticallv  stable  cr  >  <x .  4.1  [7] 

(ii)  II  G(jw)— Gr(jw)  II  ^  ^  2  Trace(CTk4.,  +  •••  +  cr, )  =  2  Trace  Z;.  [8.9]  where 

L  =  L‘  ^  .  £.  6  Rtxlc  . 


For  Hankel  norm  optimal  approximations  [&],  similar  results  hold  and  it  is  known 
that  if  G.(s)  is  a  k-dimensional  Hankel  norm  optimal  optimization,  then 

II  G(  jw)-G.(jw)  II  <  Trace  £:  (9) 

with  the  cr,  the  Hankel  singular  values  as  in  the  balanced  realization. 

From  the  above  considerations  it  is  clear  that,  while  using  the  balanced  realization  ir. 
reduced  order  modeling,  one  cannot  expect  optimal  solutions,  the  smaller  the 

c rk  +  1 . crn  the  better  the  quality  of  the  approximation  and  the  closed  one  is  to  the 

theoretical  lower  bound  II  G(jw)— Gr(jw)  II  ^  cr^]  [8]. 

FROBENIUS-HANKEL  NORM  DESIGN 

Results  relating  balanced  realizations  to  Hankel  singular  values,  and  Hankel  singular 
values  to  a  bound  on  the  H^-norm  can  be  used  to  develop  a  new  framework  for  the  solu¬ 
tion  of  many  control  problems  including  the  low-order  controller  design  problem  and  the 
controller  reduction  problem.  The  framework  is  best  described  on  the  simplest  disturbance 
rejection  problem  described  by 

x  =  Ax  +  Bu  +  Gv  ( 10) 

y  =  Cx 
z  =  Hx 

where  u  is  the  control,  v  is  the  disturbance,  y  is  the  measured  and  z  is  the  controlled  input. 
Optimal  disturbance  rejection  is  formulated  as  a  minimization  problem  involving  the  pro¬ 
duct  of  the  controllability  grammian  associated  with  the  disturbance  and  the  observability 
grammian  associated  with  the  controller  output  subject  to  admissible  feedback  controls 
from  y  to  u.  This  is  in  effect  equivalent  to  the  Frobenius  norm  defined  on  the  Hankel 
singular  values  of  the  system,  and  is  therefore  called  the  Frobenius-Hankel  norm.  The 
ensuing  approach  can  then  be  shown  to  exhibit  the  following  characteristics: 

A  Frobenius  norm  associated  with  Hankel  singular  values,  called  the  FH  norm  is  used 
to  define  the  optimal  controller. 
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Controller  order  can  be  fixed  in  advance,  resulting  in  the  design  of  controllers  of 
bounded  dimension. 

The  optimal  controllers  are  guaranteed  to  be  stable  if  there  exists  at  least  one  stabiliz¬ 
ing  controller  of  the  considered  order. 

The  optimal  controller  have  guaranteed  frequency  domain  performance  as  measured 
by  the  H^-norm. 

Computational  aspects  are  significantly  less  burdensome  than  with  all  presently  avail¬ 
able  controller  reduction  or  low-order  controller  design  techniques. 

The  approach  can  be  used  in  many  control  problems  such  as  model  reduction,  model 
reference  design,  disturbance  rejection  and  others.  The  work  was  motivated  by  research  on 
design  of  low-order  controllers  via  projective  controls  approach.  It  immediately  solves  a 
disturbance  rejection  problem  emanating  from  the  projective  controls  approach  and  can  be 
used  to  select  the  free  parameters  in  the  resulting  controllers  [3],  Design  of  decentralized 
low  order  controllers  is  a  straightforward  and  important,  extension  of  the  described  frame¬ 
work.  It  enables,  in  particular,  the  allocation  of  the  free  parameters  in  the  decentralized 
controllers  obtained  using  the  projective  controls  approach  [4].  Finally,  while  it  is  of 
interest  to  determine  the  optimal  FH  norm  solution,  most  important  in  our  new  approach  is 
the  availability  of  simple  computational  algorithm  to  find  optimal  solutions  for  constrained 
problems  emanating  from  controller  parameterization,  or  constraints  on  controller  order. 
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Decentralized  control  is  intuitively  appealing  for  vibration  control  of  large  flexible  struc¬ 
tures:  It  offers  simplified  control  system  implementations  which  basically  require  the  feed¬ 
back  of  local  measurements  to  close  the  control  loop  for  systems  that  may  have  a  few 
hundred  to  thousands  of  control  loops.  Many  decentralized  control  methods  have  been 
developed  in  the  last  decade,  typically  using  one  of  these  two  approaches.  The  first  is  an 
Information  Constraint  Approach  in  which  the  control  design  uses  the  model  of  the  large 
structure  and  seeks  a  feedback  control  solution  to  the  problem,  subject  to  constraints  on 
the  flow  of  the  measurement  information  [1,2].  The  control  action  decentralization  is  in¬ 
herent  in  any  feasible  solution  of  the  problem  satisfying  the  imposed  constraints.  In  this 
approach,  the  design  utilizes  a  model  of  the  large  structure  in  an  optimization  problem  of 
the  order  of  the  model  dimension.  Because  of  the  need  to  exercise  a  model  of  the  large 
structure  in  the  design,  the  designer  is  forced  to  reduce  the  model  order  to  cope  with 
computational  resource  limitations,  at  the  expense  of  causing  spillovers  in  the  controlled 
structure. 

The  second  is  a  Subsystem  Decomposition  Approach  in  which  the  large  structure  model  is 
first  decomposed  or  partitioned  into  an  interconnection  of  subsystem  models  3-7_.  Local 
control  designs  are  produced  by  designing  controllers  for  the  local  decoupled  subsystem 
models  which  axe  derived  from  the  large  structure  model  by  discarding  the  interconnec¬ 
tions.  These  local  designs  are  solved  independently,  thereby  reducing  the  computational 
resource  requirements  for  performing  control  designs  for  very  large  structures.  The  behav¬ 
ior  of  the  locally  controlled  structure  is  analyzed  using  aggregation  analysis  methods:  An 
aggregate  variable  is  defined  to  represent  the  subsystem’s  dynamic  interactions  with  the 
other  subsystems.  The  order  of  the  aggregate  analysis  problem  is  equal  to  the  number  of 
subsystems  which  is  determined  by  the  decomposition  scheme. 


*  Work  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by  the  Lawrence 
Livermore  National  Laboratory  under  contract  number  W-7405-ENG-48. 
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This  paper  introduces  a  new  framework  for  the  design  of  controllers  for  truss  structure  vi¬ 
bration  control  which  is  closely  related  to  that  of  the  Subsystem  Decomposition  Approach. 
The  method  developed  herein  deviates  from  conventional  control  system  design  practice  in 
which  a  dynamic  model  of  the  open  loop  plant  is  often  the  initial  data  given  to  the  control 
system  designer.  Instead,  the  controlled  plant  is  assembled  from  the  controlled  components 
in  the  control  design  process.  The  development  of  this  controlled  component  synthesis 
method  is  motivated  on  one  hand  by  the  well  developed  component  mode  synthesis  meth¬ 
ods  [8-10)  -  a  collection  of  structural  analysis  methods  which  has  been  demonstrated  to  be 
effective  for  solving  large  complex  structural  analysis  problems  for  almost  three  decades, 
and  on  the  other,  stimulated  by  the  subsystem  decomposition  viewpoint  in  large  scale 
system  theory.  Connections  bewteen  controlled  component  synthesis  and  existing  large 
scale  system  decomposition  techniques  are  established  herein  to  build  a  control  theoretic 
foundation  for  the  developed  method.  A  simple  truss  vibration  control  problem  has  been 
employed  to  illustrate  the  design  procedures,  as  well  as  demonstrating  the  potentials  of 
the  developed  method  for  controlling  very  large  dimensional  repetitive  truss  strucutres. 
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Introduction 

We  present  new  ideas  which  lead  to  feedback  control  laws  for  large  flexible  structures  which  arc 
insensitive  to  model  uncertainty.  A  pole  placement  method  is  presented  which  leads  to  near¬ 
unitary  closed  loop  eigenvectors,  and  a  new  method  is  introduced  to  design  the  control  while 
simultaneously  considering  three  competing  measures  of  optimality:  Robustness  versus  Integral 
Square  State  Error  Energy  worses  Integral  Square  Control  Energy.  We  have  established  that  the 
algorithms  are  applicable  to  at  least  moderately  high-dimensioned  systems.  In  the  present 
discussion,  controls  for  two  coupled  flexible  bodies  are  considered.  A  6x24  gain  matrix  is 
designed  to  control  a  12  modes  system  using  6  actuators.  We  have  also  developed  control  law  s 
for  the  R2P2  simulator  at  Martin  Marietta:  in  this  case  3  actuators  are  used  to  control  a  1 2th  order 
system.  Simulation  studies  indicate  that  we  have  indeed  achieved  robust  designs  without  signifi¬ 
cant  difficulties  associated  w  ith  spillover  into  the  uncontrolled  modes. 

In  this  present  discussion,  we  overview  several  key  ideas  and  numerical  results.  In  the 
references,  we  provide  details  of  the  formulation,  discussions  of  salient  features,  and  connection 
to  the  availiable  literature. 

/I  Few  Details 

Consider  a  linear  autonomous  dynmaical  system  (  having  r.  states  and  in  actuators!  with  full  state 
feedback: 

x  =  Ax  +  Bu ,  u  =  -  Gx  .  1 1 

The  closed  loop  system  is  then 

x  =  (A-BG)x  2. 

and  the  closed  loop  eigenvalues  and  eigenvectors  satisfy 

(A  -BC,  -  Xt  / )6,  =  0.  /'  =  1.2.  ...  n  (3) 
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If  we  are  solving  [he  usual  foreword  eigenvalue  problem  (given  .1,  It,  ami  (i.  find  r.\  X.  O'm.  we 
would  determine  the  eigenvalues  from  the  characteristic  equation:  detl.A  -  lid  -  A  /  1  =  0. 
However,  we  seek  to  solve  the  inverse  eigenvalue  problem  (given  A,  It,  X’s  Ce  o's.  find  (I).  we 
will  find  this  problem  easier,  in  one  sense,  since  we  will  find  it  involves  only  linear  operations. 


Notice  that  the  condition  of  Eq.  (31  can  be  re-written  in  Sylvester  form  as 

(/I  -  1)6;  =  B  ht ,  i  =  1,2, ...  n  (4) 

where  hi  ~  G  6, ,  i  =  I,  2,  ...  n  (5) 

Introducing  the  following  matrix  definitions: 

//=  [/(/  hy  . .  lin  ],  <T>=  [p;  6?  • .  On  ],  A  =  diag(?*/  .  .  .  \n  )  ((■>) 

then  Eqs.  (4)  and  (5)  become 

A  -  <I>A  =  BH  (7) 

and  //=G(D  (iS) 

Referring  to  Eqs.  (7),  (8),  consider  the  following  algorithm: 


(1)  Choose  the  closed  loop  eigenvalue  matrix  A,  and  the  parameter  matrix  II. 

(2)  Solve  Eq.  (7)  for  <I>. 

(3)  Solve  Eq.  (8)  for  6'. 

As  a  practical  matter  this  algorithm  is  oversimplified,  because  there  are  an  infinity  of  possible 
choices  for  the  matrix  II,  and  most  of  them  lead  to  poor  control  laws!  Since  it  is  difficult  to 
guess  the  matrix  II,  perhaps  the  problem  would  be  easier  to  solve  if  the  parameters  which  must 
be  guessed  had  more  obvious  physical  meaning.  This  is  precisely  the  path  pursued  in  references 
[1,  2]  where  //  is  found  to  make  the  resulting  closed  loop  modal  matrix  as  near  as  possible! in  the 
least  square  sense)  to  a  specified  target  set  of  eigenvectors.  Actually,  three  target  eigenvector 
sets  are  considered: 

(i)  The  modal  matrix  constructed!  1,  2]  from  open  loop  eigenvectors. 

(ii)  The  unitary  matrix  nearest[2]  the  modal  matrix  of  open  loop  eigenvectors. 

(iii)  The  unitary  matrix  (left  singular  vectors)  resulting!  1  ]  from  singular  value  decotn- 
sition  of  the  matrix:  [(/,(/■, ...(/  ],  where  if  is  a  basis  matrix  spanning  (.1  -/.  /)  ’//. 

The  explicit  algorithms  and  justifications  for  determining  II  corresponding  to  the  above  three 
choices  of  target  eigenvectors  are  given  in  the  references.  It  is  also  evident  that  an  initial  choice 
can  be  subjected  to  further  optimization  to  improve  some  other  performance  measure  and/or 
satisfy  prescribed  constraints.  It  is  of  significance  to  record  here  the  following  observations: 

(a)  For  cases  in  which  the  closed  loop  eigenvalues  are  to  be  only  slightly  perturbed 
from  their  open  loop  values,  then  the  choice  (i)  has  been  found  to  lead  to  the 
smallest  gains,  but  not  necessarily  the  most  robust  design. 

(b)  For  cases  in  which  some  of  the  closed  loop  eigenvalues  are  moved  deep  into  the 
left  half  plane,  then  choices  (ii)  and  (iii)  lead  to  consistently  smaller  condition 
numbers  than  (i)  for  the  closed  loop  modal  matrix,  and  usually  smaller  gain  norms 
as  well.  To  date  a  consistent  pattern  has  not  emerged  as  to  which  is  best,  both 
lead  to  robust  designs  with  small  gain  norms. 
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Multiple  Objective  Optimization 


DO 

Slate  Error  Energy:  J  s  {II)  =  j.v7  Qs x  tit  =  trace (Ps X(, ) 

<> 

OQ 

Control  Effort  Energy:  J u  (H)  =  jiirQu  u  tit  =  tracc(PHX0  ) 

o 

Stability  Robustness  Measure:  J  L,  (//)  =  *(<  I>) 
where  =  -Q1 '  (7V  0  .  / (  /vj  *  +  ?e,  ), 

p'i^-^Qji'io"  +k}) 

A  (0)  =  IOI IO  ’|  =  o(cI))/ct(cI>) 

=  condition  number:  the  ratio  of  the  maximum  and  minimum  singular 
values  of  the  elosed  loop  modal  matrix  <!>. 


^  <>  /l  (a^  } 


Multiple  Index  Optimization:  (leneration  of  Tradeoff  Surfaces 

Step  1.  Optimization  of  Primary  Objective  f  unction  ... 

Find  the  II ,  matrix  which  min  adzes  the  robustness  index  J  Oh. 


Step  2.  Evaluation  of  Secondare  Objective  function  Gradients 

Define  judicious  directions  in  the  control  yain  space  to  evaluate  the  design  whit  It  mini¬ 
mizes  the  primary  objective  function.  Establish  a  two  parameter  family  <  f  pain  vanati <  ms 
(atony  the  gradients  of  the  two  secondar  .  indices ): 

q  =  col  {111  I  •  11 2/ - hmn  )  =  col  {II).  qc  -  col  {II  e  ) 


Vs 

a  is 
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a/  v  , 
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</(«,,«„  )  =  \qe  \(  —  +  r/,v  +  a 
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Step  3.  Generation  of  the  Performance  Evaluation  Surfaces 

.4  family  of  judicious  gain  variations  is  generated  by  sweeping  (fXs.«„),  since 

(Vs  ,  Vu  )  arc  the  directions  of  steepest  variation  of  the  two  secondary  indices:  The 
surfaces  J£(ct.s ,  CLU  ),  J^(CLS ,  CLU),  J  J,CLS  .  0CU ).  and  implicitely.  JJJS  .  J  J  provide 
perspective  on  the  optimal  design  and  its  neighbors. 

Three  typical  performance  variations  are  shown  in  the  figures  below.  In  this  case,  the  closed 
loop  eigenvalue  positions  (for  a  24  th  order  system  having  6  actuators)  were  held  fixed:  so  these 
surfaces  indicate  design  freedom  beyond  pole  placement.  It  is  evident  that  even  with  fixed  pole 
locations,  significant  variations  in  all  three  indices  result  from  this  family  of  judicious  gain 
variations;  it  is  also  interesting  to  note  the  gain  region!  lower  left  quadrant)  which  is  relatively  Hat 
and  near-optimal  for  all  three  indices. 
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Robustness,  State  Energy,  tfc  Control  Energy  Variations:  Beyond  Pole  Placement 


Notes  on  Presentations  and  Discussions  in  the 
Control  Session 


Contributed  by 
John  L.  Junk ins 


Alan  Jenkins  -  Integrated  Control 

Formulations  and  algorithms  for  simultaneous  optimization  of  structures  and  controllers  were 
developed.  Continuum  and  equivalent  continuous  modeling  '  LOG”  and  "Q"  control 
parameterization.  Believes  "Q"  is  superior. 

Q.  Wouldn’t  linking  finite  element  parameters  yield  as  much  reduction  in  the  number  of 
model  parameters  as  the  equivalent  continuum  approach? 

A.  Not  sure  (several  members  of  the  audience  were  confident  that  the  answer  is  "ves".) 

Q.  Shouldn’t  the  order  of  the  compensator  be  established  early,  or.  rather  than  carrying 
through  very  high  order  designs  before  considering  order  reduction  at  the  end? 

A.  There  is  disagreement  on  this  answer. 

Q.  Do  you  really  believe  it  is  feasible  to  optimize  the  structure  and  control  system 
simultaneously?  How  do  you  deal  with  different  objectives  in  the  two  designs? 

A.  Using  multiple  objective  optimization  techniques. 


Bob  Kosut  -  Adaptive  Control 

Q.  How  many  parameters  can  you  really  adapt? 
A.  At  most,  a  half  a  dozen. 


Ed  Crawley  -  Why  use  Control  Experiments  to  learn  "what  we  don't  know" 

Recommends  performers  of  experiments  carefully  document  three  things: 

1.  What  we  expected 

2  What  we  got 

3  What  went  wrong 
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Notes  on  Presentations  and  Discussions  in  the 
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Contributed  by 
Gary  L.  Slater 


A.B.  JEN'KIM  (ISAAC) 

Q:  In  control  performance  index,  shouldn't  one  include  the  mass  of  tne 

power  supply,  etc.  in  the  overall  optimization? 

A:  Yes,  it  must  bo  included,  although  it  is  not  shown  here  and  is  difficult 

to  do . 

Q:  Analysis  uses  controller  of  fairly  high  order.  There  is  a  problem  here 

as  there  are  no  flight  qualified  hardware  of  this  type. 

A:  Mil  Spec  FIS  filters  are  'almost'  available. 

Q:  What  are  you  doing  for  structural  design?  I  do  not  understand. 

A:  °ut  areas,  masses  into  the  performance  index,  e.g.  in  a  truss  type 
structure  or  use  areas. 

Q:  How  do  you  do  the  Q-factorization?  Previous  results  show  some 
difficulty . 

Q:  Have  you  thought  about  local  minima?  How  would  you  handle  this? 

A:  Typically  start  at  different  locations  and  see  if  convergence  is  to  the 
same  solution.  Do  not  know  how  high  a  problem  this  is. 

Summary  of  a  fairly  long  discussion: 

Q:  What  is  the  correct  perf-index  for  the  control-structure  interaction? 

A:  That's  a  difficult  question. 


D.  MILLER:  Waves 


Q:  You're  talking  about  sinusoidal  waves  -  there  is  no  advantage  over  modal 

analysis.  Wave  approach  is  good  for  'discontinuous'  waves  where  you 
need  an  infinite  number  of  modes. 

A:  The  advantage  is  that  you  only  need  local  properties  nere.  Do  not  need 
modes  at  all. 

Q:  What  kind  of  sensors  are  required^ 

Q:  Implication  that  modal  analysis  does  not  work  is  incorrect! 
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G.  HYLAND:  Theory  and  Practice 


Q:  What  about  nonl inear i ties  in  joints? 

A:  Joint  nonlinearities  still  there.  Stiffness  »  l lower  it  low 
amplitudes.  Testing  shows  'modes'  do  still  exist  however. 

Q:  How  many  sensors/actuators  on  plate? 

A:  4  Actuators,  16  Sensors.  Data  not  yet  published  -  will  be  soon. 

Q:  What  kind  of  actuators? 

A:  Coil  suspended.  Rare  earth  magnet  mounted  on  _ (?) 

Q:  There  are  similar  experiments  done  on  circular  plate?  Is  square  plate 
e as i or? 

A:  'Jo  number  of  modes  is  roughly  similar.  Key  is  the  cuts  (Note:  Plate 
has  diagonal  cuts  at  corners).  Help  prevent  buckling  effects  due  to 
slight  curvature  in  plate. 

Q:  Ghow.n  that  by  changing  sensors,  you  change  the  zeros  and  therefore  the 
effect  of  the  actuators  -  Did  you  look  at  this? 

A:  We  looked  at  sensor  location  but  cannot  comment  on  how  we  got  final 
location. 


J.  YOUNG  -  Controlled  Component  Synthesis  Method 

Comment  (41)  Good  to  see  control  engineer  adapt  algorithm  to  structure  - 
rattier  than  other  way  around.  Need  to  look  at  time  response 
more  carefully.  Time  response  shown  may  have  "beats". 

(42)  This  is  addressing  one  problem  -  looking  at  behavior  of  local 
structure  -  Short  iength  scale.  Need  something  done  to  look  at 
longer  length  scales. 

(43)  This  is  a  'local'  design  technique.  There  is  a  trade  off 
between  centralized  design  and  losing  information  about  global 
structure . 


J.  JUNKINS 

Q:  Can  you  draw  conclusions  when  ft  actuators  =  4  sensors. 

A:  In  that  case  you  get  IMSC . 

Q:  3ack  to  trade-off  between  robustness  and  performance  -  What  ibout 
performance? 

A:  Fixed  gives  approximately  same  performance  (factor  of  3  or  4)  but  do 
give  great  difference  in  robustness. 
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QZGUNER  -  Decentralized/Relegated  Control 

Q:  Please  comment  on  mirrors  on  beam;  "servo- compensated"? 

A:  Small  motor  on  each  mirror;  Just  use  disturbance  model  in  feedback 
compensation. 

Q:  What  sensors  pick  up  vibration? 

A:  Use  a  tip  accelerometer. 
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Summary  of  Presentations  in  the 
Control  Session 


Contributed  by 
K.  David  Young 


!' r i »!'**> -<>r  I’mit  O/guner.  Ohio  State  \  niversitv.  outlined  the  Advanced  Beam  C  ontrol  I'.x- 
p. -riinent-  to  he  carried  out  at  AI'W'AL.  The  first  of  this  series  of  experiments  will  involve  a 
cylindrical  brain  suspended  vertically  with  a  disk  attached  to  the  free  end.  l  our  proof  mass 
actuators  are  mounted  on  the  disk  to  control  excite  the  beam  disk  structure.  Three  welj- 
lontrol  design  methodologies  will  be  selected  lor  the  validation  of  control  svMem 
di-'iun  performance  against  experimental  data  on  a  controlled  flexible  structure,  ('mit  pre- 
sented  ;i  framework  for  control  authority  relegation  -  research  supported  bv  bDIO  LLNL 
A  lielegator  is  a  key  element  in  his  relegation  scheme  -  it  serves  to  optimize  the  set  points 
for  the  subsystems  in  the  large  scale  system.  The  developed  control  authority  relegation 
scheme  has  been  applied  to  the  slewing  and  vibration  control  of  a  flexible  structure  with 
multiple  mirrors  attached,  in  particular,  to  a  slewing  mirror  structure  experiment  at  OST. 
Modeling  aspects,  such  ns  the  optical  rav  tracing,  the  integration  of  actuator  dynamics  with 
structural  dynamics,  as  well  as  unit li variable  servomechanism  design  results  were  presented. 

Professor  Juri  Medanic.  t’niversity  of  Illinois,  presented  an  overview  of  the  current  status 
oi  decentralized  control  of  large  scale  systems.  In  particular,  he  summarized  a  projec¬ 
tive  control  approach  for  reduced  order  controller  design,  and  some  results  on  a  Linear 
Quadratic  Irequency-'haped  optimal  regulator  approach,  both  of  which  were  supported  h\ 
ST'IO  I.l.NL.  These  r  wo  approaches,  when  integrated  avoid  exciting  unmodelled  dynam¬ 
ics  in  control  system  design  for  flexible  structures.  Juri  highlighted  a  new  control  design 
framework  for  disturbance  rejection  that  the  Illinois  group  is  currently  developing.  This 
framework  utilizes  a  troheniiis-Hankel  norm  to  produce  a  low  order  controller  which  allows 
the  disturbance  rejection  problem  to  be  solved  using  a  projective  control  approach. 
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